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Abstract 

Finite-size critical systems defined on a parallel plate geometry of finite extent along one single 
(z) direction with Dirichlet and Neumann boundary conditions at z = 0, L are analyzed in momen¬ 
tum space. We introduce a modified representation for the discrete eigenfunctions in a renormalized 
one-particle irreducible vertex part (IPI) scalar field-theoretic framework using either massless or 
massive fields. The appearance of multiplicities in the Feynman rules to construct diagrams due 
to this choice of representation of the basis functions is discussed along with the modified normal¬ 
ization conditions. For nonvanishing external quasi-momenta, Dirichlet and Neumann boundary 
conditions are shown to be unified within a single formalism. We examine the dimensional crossover 
regimes for these and show a correspondence with those from antiperiodic and periodic boundary 
conditions. It is demonstrated that finite-size effects for Dirichlet and Neumann boundary condi¬ 
tions do not require surface fields necessarily but are implemented nontrivially from the Feynman 
rules involving only bulk terms in the Lagrangian. As an application, the critical exponents g and v 
are evaluated at least up to two-loop level through diagrammatic means. We show that the critical 
indices are the same as those from the hulk (infinite) system irrespective of the boundary conditions. 

PACS numbers: 64.60.an; 64.60.F-; 75.40.Cx 


* e-mail:messiasvilbert@df.ufpe.br 
^ e-mail:jborba@petrobras.com.br 
^ e-mail:mleite@insti.physics.sunysb.edu 


1 



I. INTRODUCTION 


Fields confined in a certain space region have especial properties whose study is worth¬ 
while. Focusing entirely on conventional critical behavior, they correspond to the order 
parameters of certain systems undergoing a second order phase transition [i|, 1^. The or¬ 
der parameter has well dehned values at the boundary surfaces. This can be implemented 
through boundary conditions, e. g., periodic (PBC), antiperiodic {ABC), Dirichlet {DBC), 
Neumann {NBC), etc. and/or via external surface fields. 

A layered parallel (hyper)plate geometry is dehned by {d — l)-dimensional (hyper)planes 
of inhnite extent whose limiting plates are located, say, at z = 0,L {z is the space direction 
perpendicular to the “planes”). We can either vary the distance between the boundaries or 
keep L hxed and study universal amounts for the conhned system. The change of funda¬ 
mental quantities, like energy, with the variation on L originates the Casimir effect 3l-l6|, in 
which long range forces arise in response to the conhnement of huctuations. A similar effect 
in conhned geometries of this type is the thinning of ^He near the superhuid transition due 
to a force causing the change of the surface free energy in the critical huctuations whenever 
the order parameter vanishes {DBC) on both boundary surfaces [^. 

By maintaining L hxed, however, the critical behavior of the order parameter in thin hlms 
could manifest hnite-size corrections. BBC resembles closely experimental conhned systems 
as in hows from glassy materials above and below the jamming transition in comparison with 
the unconhned behavior, whereas DBC for the same system indicates a surface phenomenon 
js]. Nevertheless, DBC and BBC take into account purely hnite-size contributions in the 

between parallel 


experimental discussion of conhned ‘^He 


9|] as well as other huids jl^ . 
plates. These corrections also take place in the measurements of spatial correlation of current 
critical points in an open billiards system [n| in describing the helds of turbulence, sound 
waves and acoustics, among others. On the other hand, Neumann boundary conditions 
govern the transition bulk-surface when the critical behavior deviates from the {N, d) to the 
{N,d — 1) universality class and a nontrivial mixing of hnite-size and surface ehects was 
argued to take place [l2|. Moreover, DBC and NBC have been investigated on a wide class 
of hnite-size models, like amorphous solid suhering plastic deformation in a certain class of 
depinning models 13|, dynamics of order reconstruction in conhned nematic liquid crystals 


14j |. determination of free energy corrections in the conhned spherical model 1^, etc. 
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Phenomenological finite-size scaling theory states that close to the bnlk critical temper¬ 
ature, the variable j (£ is the bulk correlation length) measures the deviation from the 
bulk critical behavior 


16 


20|. According to some authors (l6l-ll8l|. the description of the 
finite system is not limited to the values ^ they suggested that the finite-size critical 
exponents should be identical as those from the bulk (infinite) system. 

The description of finite-size critical systems using momentum space renormalization 
group field-theoretic methods has basic aspects which are simple to grasp. First, for parallel 
plate geometries the boundary conditions on the plates are implemented through the bare 
free propagator. Second, the typical length L separating the boundaries can be included 
in the Feynman rules. For PBC and ABC, the momentum along the finite-size direction 
turns into the quasi-momentum and the integral along the 2 : direction gets transformed to an 
infinite sum (Nemirovsky-Freed {NF) method 2l|]). Motivated by the NF Green’s functions 
formalism, a one-particle irreducible {IPI) vertex part framework was designed recently 
to the treatment of finite-size systems subject to PBC and ABC using either massive or 
massless fields. Within this finite-size technique, scaling theory holds rigorously in the whole 
region 0 < j < 00 in agreement with previous claims 18|, l20|. Explicit computations at 


higher order were performed corrobora 


ing that the finiteness corrections are not sufficient to 


modify the bulk critical exponents 22(]. A different piece of folklore exists in the literature 


regarding DBC and NBC: they represent free surfaces and are appropriate to describe 
finite-size plus surface effects. Indeed one can introduce surface fields in conjumination 
with bulk fields for this boundary conditions breaking translation invariance in this way 
23l |. If we allow only bulk fields, ruling out ordering surface fields associated with boundary 
plates, is that possible to renormalize solely bulk fields subject to those boundary conditions 
which do not produce any surface contributions, but with manifest translation invariance 
breaking? Can a framework for DBC and NBC be devised to explore the full finite-size 
region 0 < j < ool If we achieve consistency with the results obtained from PBC and ABC, 
the emerging description represents a modern finite-size scaling regime. Is there any simple 
relation with the associated unconfined system (the “bulk” criticality)? How the decreasing 
of the confinement region rules over those situations? How the presence or absence of surface 
fields alters the criticality? 

In this work we build up an ab initio renormalized one-particle irreducible vertex part 
formalism in momentum space for Dirichlet and Neumann boundary conditions in order 
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to calculate perturbatively universal critical properties including only bulk fields. For sake 
of comparison with previous one-loop results, we first employ massive fields. We modify 
the Feynman rules with respect to the previous approach using sines and cosines as basis 


functions 


2l| due to our decomposition of them in terms of exponentials only. The Feynman 


diagrams obtained from those rules can be expressed in terms of integrals (and summations) 
identical to the PBC and ABC cases, but contain extra “nondiagonal” terms in which there 
are not as many summations as there are integrals over the {d — l)-dimensional transverse 
space. 

This new feature will permit us to prove that each diagram is composed of two parts: 
the first one, where the momenta as well as quasimomenta are conserved in every diagram 
order by order in perturbation theory and the second, substantiated by the presence of the 
“nondiagonal” terms which break the translation invariance (violation of quasi-momentum 
conservation) along the finite-size direction. These terms are important for they make a 
clear-cut distinction in comparison with the simplicity obtained from PBC and ABC re¬ 
sults. They look like surface contributions at first sight albeit they are purely finite-size 
corrections and perhaps this is the mathematical origin of the folklore above mentioned. We 
will show that they do not contribute to the leading singularities in dimensional regulariza¬ 
tion. We introduce new Feynman rules for vertex parts including composite operators which 
can be renormalized multiplicatively. We demonstrate that for a certain quasi-momentum 
distribution of external legs not belonging to composite fields, the external quasi-momentum 
of the insertion of the composite operator admits more than one combination of the external 
quasi-momenta of the other usual external legs (not associated to composite operators) and 
should be properly taking into account in the set of rules order by order in perturbation 
theory. 

We discuss the dimensional crossover criteria and show the consistency with the previous 
results from PBC and ABC in 2^. Staying away from the problematic region where the 
e-expansion ceases to give meaningful results, all loop integrals considered will be shown to 
have the general structure made out of bulk plus finite-size terms with the latter depending 
on the boundary conditions. Within the finite-size plus bulk regime, we choose to renormal¬ 
ize the field theory with zero external momenta and nonvanishing external quasimomenta 
in order to unify Neumann and Dirichlet boundary conditions in a single framework. After¬ 
wards, the unifying formalism of massless fields using nonvanishing external momenta and 
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quasi-momenta for DBC and NBC is presented along with a discussion of the dimensional 
crossover criterion where the hnite-size regime starts to give meaningless results. We prove 
its equivalence with the massive case. As an application we compute the critical exponents 
rj and u using diagrammatic means. We show that the universal results are independent of 
the boundary condition. 

In Sec. II we discuss how the discrete eigenfunctions corresponding to NBC and DBC 
are expressed in terms of exponentials and the consequent modihcation of the Feynman rules. 
We construct the tensor couplings for all the primitively divergent vertex parts which can be 
renormalized multiplicatively and demonstrate the emergence of multiplicities with respect 
to the bulk (inhnite) theory due to the exponential representation of the basis functions. We 
show that DBC possess diagram multiplicities identical to NBC whenever both theories are 
dehned at nonvanishing external quasi-momenta. In this picture we compute explicitly some 
sample Feynman diagrams. We conclude by presenting the unified set of all diagrams which 
are going to be relevant in our discussion. The arguments are valid for both formulations of 
massive and massless helds. 

We discuss the renormalization of massive helds in Sec. III. We start with a large number 
of diagrams and achieve the reduction to a smaller set of graphs owing to the nontrivial 
cancellation of the mass insertions. This feature includes both the diagonal and nondiagonal 
terms of tadpole diagrams. This argument is decisive in proving that the theory can be 
renormalized without making any reference to ad hoc surface helds. Normalization conditions 
are dehned for the primitively divergent vertex parts in order to assure the hniteness of all 
vertex parts that can be renormalized multiplicatively. We briehy discuss how the how in 
the mass scale ahects the renormalized vertex parts in essentially the same way as in the 
bulk case. The limit L —)■ cxd is shown to retrieve the bulk result, whereas the L —)■ 0 limit 
marks the onset of the dimensional crossover which invalidates e-expansion results. We also 
list the solution for the higher loop diagrams from Appendix A which will be required in 
the determination of critical exponents in Sec. IV. The dependence of the renormalization 
constants/functions on the boundary conditions disappears in the hnal expression for the 
exponents. 

Section V contains the explicit discussion of massless helds and their multiplicative renor¬ 
malization for both boundary conditions in the unihed description. The examination of how 
mass insertions are cancelled is explained in two different ways in order to get a minimal 
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number of diagrams to work with in the determination of critical exponents. An in-depth 
discussion on the validity of the hnite-size regime with the approach to the bulk criticality 
as well as to the dimensional crossover regime is presented. We give a brief description of 
the solution to the higher-order massless diagrams by writing down their expressions and 
point out their similarity with PBC and ABC arguments for the massless helds. In Sec. 
VI we compute the critical exponents in the massless approach. 

In Sec. VII we discuss our results. The dimensional crossover regime is analyzed explicitly 
by focusing on the one-loop correction to the bulk case of the four-point vertex part. We 
perform a comparison of these regimes with those from PBC and ABC previously studied. 
We vary the correlation length and establish that even in the massive theory there are regions 
for hxed, hnite ^ where j < 1 and the e-expansion results are still valid. 

Section VIII displays the conclusions and future directions within the formalism intro¬ 
duced in the present work. In addition, we point out how these new aspects can be adapted 
to tackle the problem of competing systems. 

In Appendix A we compute higher order massive integrals. We decided not to give a 
detailed account of the solution of massless integrals of this problem in another appendix 
for the same reason. The reader is advised to consult Appendix A and Ref. for grasping 
the details. 


II. MODIFIED NF APPROACH TO FEYNMAN RULES WITH DBC AND NBC 


A. Review of NF approach 


First we will discuss briefly the held theory setting introduced by Nemirovsky and Freed 
{NF) 2l|] for constructing Feynman diagrams in momentum space with the boundary con¬ 
ditions of interest in the present work. We are going to restrict ourselves to the situation 
where no external surface helds are allowed. 

The bare Lagrangian (free energy) density is composed by scalar helds with 0{N) sym¬ 
metry, dehned on the volume enclosed by the two limiting {d — l)-dimensional parallel 
hyperplates located at z = 0, L (bulk helds). It is given by: 


— 2I V 00^ + 


( 1 ) 


where (po, /iq and Aq are the bare order parameter, mass (/Xq = to is the bare reduced 
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The 


temperature proportional to '^^r) coupling constant, respectively { 2 , 

0{N) symmetry of the bare order parameter means that it is a vector of N components 
((^ 0 )^ = + + The additional indices related to the 0{N) symmetry of the fields 

were not written since they are identical to the infinite system. The space directions split 
in two distinct sets denoted by p representing the coordinates along the {d— l)-dimensional 
subspace parallel to the plates and the perpendicular 2 ; axis, collectively represented in the 
vector form as a: = {p,z). The space is layered and filled with parallel plates in the region 
between z = 0 and z = L. The held satishes (j)o{z = 0) = (j)o{z = L) = 0 for Dirichlet 
boundary conditions (DBC), whereas ^^{z = 0) = ^^{z = L) for Neumann boundary 
conditions {NBC). 

The order parameter can be related to its Fourier modes in momentum space through 

4 >q{x) = I d^~^kexp{ik.f>)Uj{z)(f)Qj{k), where k is the momentum vector characterizing 
j 

the {d — l)-dimensional space. The basis functions Uj{z) satisfy the differential equation 

— = KjUj{z), Kj being the quasi-momentum along the z-direction. The eigenfunctions 

are orthonormalized according to '^Uj{z)uj{z') = 6 {z — z') and dzuj{z)u*,{z) = 6 jji. 

j 

Here Kj = ^ = aj, where j = 1,2,... for DBC and j = 0,1,2,..., for NBC {d = ^). 
We can attach a label for each boundary condition such that r = —1 corresponds to DBC 
and r = 1 to NBC. The free bare massive propagator (/Xq 7 ^ 0) in momentum space for 
either boundary condition is given by the expression cl^j{k,j) = Although not 

explicitly mentioned so far, it is obvious that the eigenfunctions depend upon the boundary 
conditions. We will suppress their relation with the boundary conditions but can retrieve 
this dependence whenever it is appropriate. 

A generic Feynman diagram involves the product of many bare propagators G^^j con¬ 
tracted with interaction vertices. The amazing fact about the structure of the hnite- 
size to take effect on the field theory is that it is implemented as an internal symme¬ 
try. For example, beside the standard tensorial couplings of the inhnite theory corre¬ 
sponding to the A^-component order parameter, each momentum line (propagator) must 
be multiplied by dzuji{z)uj2{z) and the vertices are multipled by the tensor 

^ 31323334 . ~ lo dzuj^{z)uj 2 {z)uj^{z)uj^{z). Forgetting about the composite operators for the 
time being, this is nothing but the enhanced internal symmetry representation of the di¬ 
rect product 0{N) X (/s), where the (/s) symmetry is represented by the tensors and 
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Sj^j 2 j 3 j 4 ,- The basis functions can be written as ^^( 2 ;)= ^sm(Kj 2 ;) for DSC, whereas 

1 

for NBC we have the nonzero mode as in DBC u^^~^\z) = cos{Kjz) (j = 1,2...) as 

well as The Feynman rules for vertex and propagators are given in 

In the computation of Feynman diagrams as the set of integrals in (d — 1) dimensions 
in conjunction with inhnite summations, it is difficult to establish a direct comparison with 
periodic (BBC) and antiperiodic boundary conditions {ABC). The reason is that j > 0 for 
DBC and NBC, but varies in the interval (—cxo, cxo) por BBC and ABC. The simplihcation 
which took place for BBC and ABC when we could compute integrals with all external 
quasi-momentum set to zero no longer occurs for DBC. In any serious attempt to unify 
the framework for DBC and BBC, we should hgure out how to compute graphs with 
nonvanishing external quasi-momentum for those conditions. 

Fortunately, at an external quasi-momentum symmetric point, we can compute the di¬ 
agrams provided some modihcations are introduced in the moding of the label j and ad¬ 
ditional trivial orthonormality properties. Next, we shall introduce new notation in order 
to give a unihed description of the Dirichlet and Neumann problems for nonzero external 
quasi-momentum. 



B. Exponential representation and unification of the Feynman rules for nonvan¬ 
ishing external quasi-momentum 


1 . One-loop diagrams for the two- and four-point vertex functions 


The construction of Feynman diagrams for DBC (sine) and NBC (cosine) takes into 
account solely the internal structure provided by the tensors 0 {N) 

underlying symmetry appears in exactly the same form as in the inhnite L limit and we 
shall simply attach the symmetry factor to each diagram under consideration. However, 
we are going to discuss explicitly how the combination of products of the hnite size tensors 
shows up in particular one-loop graphs. We wish to express them in terms of summations 
varying in the range (— 00 , 00 ) and {d — l)-dimensional momentum integrals. 

Let us write the basis functions as 


/ 2 \ 2 

r 1 , ... J 



\L) 

_2i ^ _ 
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{DBC), 


( 2 ) 





{NBC). 


(3) 






We can compute the tensors in this new representation. For instance, using the dehnition 
of and restricting ourselves only to positive values (ii, j 2 ^ ZC we achieve the following 


3132 

unihed form 


cF) _ 


- 1 ) 




2L 


_j_ g-*5-(ii+i2)2 _|_ ^^ia{ji-j2)z _|_ ^g-i5-(ji-j2)2j (^4^ 


which by the change z —)■ — 2 : in the second and fourth terms turns out to be given by: 


~ ^ ['''^ii-i2,0 + ^ji+j2fi] = (-1) ^ - 32 ) + h(ji + j 2 )]- (5) 


The new notation for the Kronecker’s delta is going to be very useful in what follows. We 
have to be careful with the moding since NBC allows ji,j 2 = 0 in the construction of a 
generic Feynman diagram: in the contraction of these tensors, this value of j can occur 
and all components of the hnite-size tensors should be computed (as the internal quasi¬ 
momentum indices are summed). For all practical purposes in the present paper, the result 
for the two-index tensor is simply = S(ji — j 2 ). 

It is a simple task to evaluate the tensor with the arguments at hand. We hnd 


cW 


^{ji + 32 + js + 3^) + ~ 32 + is + J 4 ) + 'rh(ji -|- j 2 — is + ji) 

+ + j2 + is ~ i4) + h(ii — i2 — is + i4) + h(ii — j2 + is — i4) 

+ <^(ii + 32 - is - i4) + r6{ji - is - is - i4) • 


( 6 ) 


In PBC and ABC the quasi-momentum conservation is represented by the condition 
^jij 2 j 3 j 4 ~ -^ 4 )’ whereas DBC and NBC in the above expression possess 

the quasi-momentum conservation as well as additional contributions as expected from the 
intrinsic nontrivial nature of those boundary conditions. 

The last expression can be further simplihed utilizing the notation j± = iiiis, j± = i 3 ±i 4 
which yields 

-^£1^4 = ^ (7) 

Q:i=±l 

Q:2=ibl 

Let us consider the components of the hnite-size tensors when j = 0. In the simplest case 
of the two-index tensor, it is easy to show that zero modes do not mix with nonzero modes. 
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which can be expressed in a compact form as When we analyze the fonr-index 

tensor, the only difference with respect to the components with all nonvanishing snbscripts is 
different normalization factors. It is straightforward to prove that < 5 'oj]j 2 i 3 ~ “*“-^2 + 

j3)+S{jl-j2+j3) + S{jl+j2-j3) + S{jl-j2-j3W\ ^2^2 = Z [<^0'!+^ 2 )and 
S'oooji = If is important to mention that neither nor are eqnivalent 

to NF notation. The above tensor is similar bnt has different components when compared 
with the NF connterpart (obtained in terms of sines and cosines) as can be verified in Eq. 
(A 6 ) from the second paper in Ref. Therefore, the components of the fonr-index tensor 
in the exponential representation enlarges the possibilities for the finite-size indices. 

These elements will snffice to onr constrnction of vertex parts diagrams not inclnding 
composite helds yet. Before working out explicitly the one-loop graphs, let us introduce 
the notation = 27rS',-(^- ■ • and from now on we are going to construct the integrals 

associated with Feynman diagrams with this modified tensor. Each loop integral will be 
represented by a (d — l)-dimensional integral multiplied by summations involving products 
of and This has the virtue of producing the metric factor a multiplying each 

loop integral. 

In practice, we have to use the tensor structure to build up the diagrams perturbatively. 
The coupling constant factors will be omitted in all diagrams to be discussed, but will be 
retrieved during the discussion of the diagrammatic expansion of each primitively divergent 
vertex part. Within a IPI vertex part framework, those primitively divergent are the 
two-, four- and two-point with one composite operator vertex parts, which are represented 
by T*^^) and respectively. Whenever we do not mention explicitly the external 

quasi-momentum index in a certain graph, these modes are arbitrary. For instance, a generic 
graph from T*^^) has external quasi-momenta labeled by R and Z 2 , T*^^) has associated external 
quasi-momenta ^ 1 ,^ 2 , ^3 and and so on. On the other hand, we are going to attach them 
to the particular diagram whenever we choose a particular fixed distribution of external 
quasi-momenta modes. 

As a hrst sample computation, consider the one-loop “tadpole” contribution to T^^^. Its 
diagrammatic expression reads 


a 


(iV + 2) 


E qir) q 


(" t ) air) 
W2hl2 


d'^-^qGoiqJi). 


( 8 ) 




Let us consider explicitly Neumann boundary conditions (r = 1). The summation splits in 
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the contributions of zero and nonzero modes as: 


OO n n OO 

lJ 2>0 ^ ^ 1=1 


,^ 2^0 


d-‘-'kG„(q,l)+ Y. S'ild.A I d'‘-'kG„(q,h) 


( 9 ) 


h,h=^ 


Note that Go{q, 1) = Go{q, —1). The term proportional to vanishes. Replacing explicitly 
the values of the tensors already derived and dehning i± = ii± Z 2 , we hnd: 


W gi-r) I ^d-li 




+6{i+ - 21) + 25(4) + d{i- + 21) + 5{i_ - 21) + 25(z_))G'o(g, /) 


Next, the identities 


2(5(4) + d(i-))Go(q,0) + ^(5(4 + 2/) 


1=1 


( 10 ) 


5(^± + 20Go(g,0= ^d(z± + 20Go(q,0-d(z±)Go(q,0), (11a) 

1 ^ 0 , l=—oo l=—oo 

OO OO 

d{ii)G{q, 0) + 25(z±)5^Go(g, /) = 5(z±) Go{q, /), (11b) 

/ = 1 l = — 00 

will be useful in order to achieve our goal of transforming summations with a bounded lower 
limit into those with unlimited negative values for the index 1. After some manipulations 
the hnal expression can be rewritten as 


a 


= a- 


{N + 2) 


» OO 

J ^ 


[5(i_) + 5(4) T <^(4 T 21) + 5(z_ — 2/)] 


l=—oo 


4 + 5-2/2 ^2 


( 12 ) 


For DBG boundary conditions we have instead 


a 


= a- 


iN + 2) 




(13) 


^ 1,^2 = ! 


With li,l 2 > 0, = 5(/i — 12 ). We can perform one summation, say over I 2 and perform 

the change of index /i —)■ / in the remaining summation. Focusing only on the summation 
and replacing the expression for we get to 

CXD „ CO 

/ d^-^kGo{k, 1) = d^[-25{i+) + 5(4 + 21) + 5(4 - 21) + 25(G) - 5{t_ + 21) 


1=1 


1=1 


-5(G -2/)] / (f-^kGo{k,l) 


(14) 
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From the identity ^((5(z±+2/)+5(z±—2/))Go(?! 0 = S 5(i±+2/)Go(q',/) S'lid by employing 

^ = 1 / = —OD 

Eqs. CD. we obtain 


a 


= a- 


iN + 2) 


/ oo 

<^“-'1 E 

1 _ 


/ = —OO 


[5(z_) - (5(4) + (5(4 + 21) - 6{i_ - 21)] 
+ 0 - 2/2 ^2 


(15) 


Comparing the DBC (r = —1) and NBC (r = 1) results, we can write down the expression 
corresponding to this graph in a unified fashion for these boundary conditions as 


a 


= a 


{N + 2) 


» oo 


l=—oo 


[(5(z_) + r(5(4) + (5(4 + 2/) + r(5(i_ - 2/)] 

4 + 5 - 2/2 ^2 


(16) 


The summations in the third and fourth term can be performed and we are left with 


a 


(A^ + 2) 


((5(f_) + r(5(4))T>i(cT,/io) + rT)i(y,(T,/io) + T»i(y,cx,^o) ,(17) 


where 

~ f d!^~^ Q 

= (18b) 

Although the first and second terms in the diagram expression have a clear counterpart in 
BBC and ABC, the third and fourth terms are entirely new, since no summation appears 
there. They look like surface terms albeit coming purely from finite-size contributions: 
they reflect the breaking of translation invariance for DBC and NBC along the finite-size 
direction. In what follows diagrams of this type (“tadpoles”) will not be required in the 
present method as we are going to show in the discussion of normalization conditions. 

Consider the one-loop contribution of the four-point vertex part T^'^b Its graph in the 
present setting is given by 

) 0 ( = 14^ £ SXUlNC. f ‘‘“-'qG.U + k, h)Go(q, k), (19) 

where k is the external momentum, ii,i 2 ,i 3 and u are the external indices from quasi¬ 
momenta. We start with r = 1 [NBC). Focusing only on the summmation, integrals and 
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propagator products, we first decompose the whole thing in the following form 


E I d‘-\,G„(q + k,k)G„(q,h) = j rf-'V;[S'*'+ fc, 0) 

lul2>0 ^ ^ 

oo oo 

Go{q, 0) + Gq{ q + k, li)Go{q, 0) + S^lxih^ohLuGo{q + k, 0)G'o(g, k] 

ll = l /2 = 1 


c(i) 


- E Sl:LM^s^M^ + k,h)Goiq:l2) 

ll,l2 = i 


( 20 ) 


Note that the single summations are actually equal. For the sake of simplicity, the dehnitions 
i± = ii ± *2 and i± = will have their utility in what follows. Using the expressions 

for the tensors and the symmetry of the integrals yield 

OO n n 

E / d'‘-'qGo(q + k,l,)G„(q,k)=4a^ / d‘*-*,(|i(4) + [<S(4) 

OO 

+(5(L)]G'o(g + k, 0)Go{q, 0) + E['^(*+ + 0 + <^(*+ - 0 + S{i- + 1) + - l)Mi+ 

1=1 

^ OO 

+/) + S{i+ -l) + h(L + /) + 6{i_ - l)]Go{q + k, l)Go{q, 0) + ^ E + ^ 2 ) 

li,h=l 

+(i(z-)- — /i + I2) + <^(*+ + — ^2) + — ^2) + 5 (f— + + I2) + h(i_ — /i + I2) 

+5(f_ + /i — I 2 ) + 5(f— — li — l2)][S{i-\- + /i + I 2 ) + + I 2 ) + + /i — I 2 ) 

+( 5 (f-)- — li — I2) + S(^i— + /i + I2) + 6 (^i— — + I2) + 5 (f— + /i — I2) + h(i_ — /i — I2)] X 

Go{q + k,li)Go{q,l2)^ ■ (21) 


The summation involving the combination [(5(f± + 1) + 5(f± — /)] can be transformed into one 
whose range extends to negative integer values not including the zero. Indeed, the identity 

OO OO 

E[^(^± + 0 + Gg - /)][h(*+ + /) + S{G -l) + S{i- + l) + h(L - /)] = E 

1=1 l^0,l=—oo 

+1) [h(?+ + /) + 6{i+ -l) + 6{i. + l) + h(L - /)], (22) 


can be used so that 


/ 0 - l)MG + l) + Gg - 0 + + 0 + - l)]Go{q + k, 1) 

^ 1=1 

/ OO 

d‘'~"q E + ^) + + ^) + + ^) + 

l=—oo 

Go(5'+ A:,/)Go( 5', 0) — 2[h(i+) + (5(f_)] [5(f+) + h(i_)] f df ^q'Go(q'+ A:, 0)Go(q', 0). (23) 
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Using these expressions inside the double summation and carrying out the computations for 
it in a similar way, we obtain the result in terms of summations with indices varying in the 
range (—cxo, cxo). The single summation is cancelled by the zero mode of each index in the 
double summation and we are left with the following expression for this diagram 

+ h) + + / 2 )][<^(*+ + + h) + — li 

Il^l2 = — O0 

+^ 2 ) + ~ ^ 2 ) + ^{'i+ ~ h ~ ^ 2 ) + + /i + I 2 ) + h(i_ — /i + I 2 ) + S{i- + /i — I 2 ) 

+5{i-— li — I 2 )] j + ^) ^i)G'o(9) ^ 2 )- (24) 


This form can be further reduced in such a way that at maximum one summation for each 
loop graph takes place. We can perform the summations very easily and express the diagram 
in the simpler form 


o 


(A^ + 8 ) 
9 


a< [d(4 + 4) + (5(4 - 4 ) + (5(4 + i_) + (5(4 - i-)]l 2 ik, 4; a, po) 


+ [(5(z_ + 4) + d(i- - 4) + d(i- + i-) + S{i- - i_)]/ 2 (/c, i-] a, po) + hik, 


4 T 4 *+ — *+ 


~ \ Til ^“t“ ^“1“ ^“1“ ~ \ r / 7 ^^— ~ \ r / 7 

^) ho) + -' 2(^7 —^ 1 ho) + -' 2(^5 —-—, —-—, ( 7 , fio) + hyk, 


2 ’ 2 
4 — i- 


---, a, po) + h[k, ---,---, cr, po) + h[k, ---,---, (T, /iq) + hik, 


i_+f_ „ , -r ,, - Z_ Z_ + i_ „ . 

7 7 ^7 ho) T 4(^7 7 7 *^7 ho) 


( 25 ) 


where 


00 « 

/2(A;,b5-,^o) = d ^ / 

/=-oo 

h{k,i,j,a,iJ,o) = d 




[g 2 + d 2/2 + ^l][(q + h )2 + + *)2 + ^, 2 ] > 


[g2 + 0-24 + p^][(g + A:)2 + + pg] 


(26a) 

(26b) 


For DBC, the absence of zero modes in the four-index tensor makes the computation easier. 
Now the diagram reads: 



(A^ + 8) 

9 


X 2 




(-1) 5(-l) 


d^ ^gGo(g + ^7 /i)Go (97 ^2)7 


(27) 


and we have to transform the summation to the range (—cxo, cx)) just as before after using 
the tensor components for r = —1. The unihed form holding for both r = —1 {DBC) and 
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r = 1 {NBC) can be written as 


o 


(A^ + 8) 
9 


a< [5(4 + 4) + 5(4 - 4) + r(5(4 + i_) + 5(4 - i_))]l 2 {k, 4; a, /xq) 




+ [r(5(z_ + 4) + <5(i- - 4)) + 5(i_ + i-) + S{i- - i_)]/ 2 (fc, a, /io) + r ^{k, - 

—^/^o) + h{k, —^^—, CT, /io) + /2 (k, —^^—, a, /io) + Hk, 

i-- i- i- + i- „ .] ~ i_ + 4 *-“ *+ ~ N ?/,*-“ H *- + *+ ~ N 

- ) cr, /io) + hik, , , cr, /io) + hik, , , cr , /io) 


'^h{k-, -r-,-r-, <j , /io) + h{k-, -r-,-r-, <J , /io) 


(28) 


with the integrals dehned by Eqs. (]26|1 . From this simple analysis, we conclude that the 
Feynman diagrams in the exponential representation involve more complicated objects like 
the “nondiagonal” integrals J 2 (A:, i, d,/io). Consequently, each momentum loop integral in 
the hnite system for DBC and NBC cannot be obtained from the inhnite system through 

OD 

the substitution J d‘^k —)■ cr f d!^~^k as in BBC and ABC. The higher the loop graph, 

j=-oo 

the lengthier is the computation of the contractions of ,■ giving the particular diagram. 

J±J2JoJ^ 

Nevertheless, the same structure is preserved as in the simple one-loop examples just worked 
out. We postpone the presentation of all relevant higher loop diagrams after our next task, 
namely, the discussion of vertex parts that can be renormalized multiplicatively and include 
composite operators. 


2. The vertex tree-level and one-loop graph 

The most nontrivial feature of the construction proposed in the present work manifests 
itself in the diagrams of the composed field. The reason for that are the basis functions of 
the composite helds associated to the vertex part (and its descendents that can be 

renormalized multiplicatively), which consists of one single type for Dirichlet and Neumann 
boundary conditions. 

Applying the operator (y^ — /^o) the composite held {(f)^Q \x))‘^ and using the held 
equations (y^ — — 0, we hnd 

(V" - /^o)(4"^(a^))" = 2(y4"Ha^))"- (29) 
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The composite field also satisfies DBC and NBC, respectively, namely 

= 0 , 

^{(t'o^^\x)f\z=o,L = ^- (30) 

We decompose the composite field in terms of its components in momentum space as 

((;/)g’'^(a;))^ = f df‘~^pexp{i'p.p)v!'p{p, j)Y. Our goal is the determination of the 
j 

new basis functions of the composite field u^p{z) preferably in terms of the previous basis 
functions expressed in terms of sines and cosines. In order to achieve that, we first work 
out the representation of the composite field in connection with the product of two single 
fields computed at the same point, apply the differential operators over them and examine 
the consequences. Afterwards, we compare the results of the same operation performed on 
the above definition of the composite fields and the basis functions u^p{z). 

Using the representation of the field and taking the product of two fields at the same point, 
we find (f)o{x)4>o{x) = I kiexp(i(ki + k 2 ).p)Uj,^(z)uj 2 (z)</)oji(ki)(^oj 2 (k 2 )- Recalling 

J1 J2 

that ^ applying this operator to the product of two fields and inserting this 

in above equation, we find 

r o du^l} . , 


^ - {pI+pD - 2ho 


dz dz 


Restricting ourselves only to nonvanishing values for ji, j 2 , we can show that {z)u"^_^{z) = 

~2 ■ ■ 

x[cos{a{ji-j 2 )z)+Tcos{a{ji+j 2 )z)]. This implies that ^ = ^-^[cos{a{ji- i 2 )z) - 
Tcos{a{ji + j 2 )z)\. Using the definition Up^j 2 = {z)uj_^{z) into the above equation leads 

us to the identity {p\ + pi + + d‘^j‘1 + = 0. These manipulations suggest the 

choice for the basis functions of the composite field: 

u^j{z) = ^cos{djz), (32) 

which implies that We can now determine the tensor associated 

to the vertex part . The typical object responsible for the appearance of these diagrams 


corresponds to 


Ip / dz(pi[\x)(j)^Q\x){(j)^o\x)f. 


Replacing the Fourier expansion of the fields and composite operator, we obtain 


_p(2,l) _ 


ZbA/ ““-^1 


pid'^ W* V"* \pi+P2+p}(f)opipi)(j)oj2iP2){(j)ojil^f, (34) 
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where 


ch) 

‘^hhj 



( 35 ) 


^0 

The constant A can be determined by imposing the normalization condition 
Uj(z)uj/(z)dz = — f) + 6{j + f)] with j,f ^ 0, which implies A = L. There¬ 

fore, it can be checked that 


“ A + 32 ) + d{i + ji - is) + rh(j + ji + is) + T6{j - ii - is)]. (36) 


We can invert the Fourier transform from the product of two helds computed at the same 
point, using the identity First, we multiply it by and 

integrate over After that, multiply the resulting expression by Uji{z) and integrate 

over z using the orthonormality conditions of Uj{z). We obtain: 

(po\pJ) = ^ J d'^~^p[(po\x)]‘^exp[-ip.^ Uj{z)dz. (37) 

Replacing this result in the left-hand side of the expression which dehnes the Fourier trans¬ 
form of the composite held (with only one summation), we get to the closure relation 
j'^jUj{z)uj{z') = 6(z — z'). This completes the basic properties of the basis functions 
of the composite held. 

Before analyzing the situation for NBC let us make an important remark concerning 
the construction of arbitrary loop diagrams for the composite held. First, observe that even 
though this tensor is similar to the four-point vertex part tensor the similarity is 

not complete. We would like to make a connection between the tensors • and , 

or in other words, the graphs of with those from F*^^) by identifying, for instance two 

external legs of the latter with the composite held insertion, as is usual for inhnite systems. 

The fact of the matter is that a simple way to link the two types of graphs which gen¬ 
eralizes the situation of bulk systems is to perform identihcations between the external 
quasi-momenta. If we set = f, we hnd = d[6{2f -F ji -F ja) + 2r5(ji -h ^ 2 ) + 

S{2j' - ji - j 2 ) + T5{2j' - ji + j 2 ) + TS{2j' + ji - j 2 ) + 25{ji - ^ 2 )] (no summation in the 
repeated index). 

Looking at Eq. (]36|1 . it is tempting to perform the identihcation 2j' = j in that formula, 
but this does not solve the problem completely, for the additional contributions showing up 
in (proportional to 5(ji±J'2)) are absent in Since the moding of the composite 

operator index can take zero values for both NBC and DBC, the construction of the 
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graphs with the tensor needs to be modified. In fact, the right combination to form an 
arbitrary loop graph shonld inclnde the j = 0 and j 0 pieces components from To be 

precise, an arbitrary graph always contains the contribntions proportional to 
{j = ±(ji +^ 2 )) contracted with tensor stemming from the pertnrbative expansion 

in the conpling constant. In terms of the diagram the combination that shonld appear is 



I(i=±(ii+i2)) T^ „l(i=o)’ 


(38) 


where the sign of j is fixed. With this recipe ntilized in the constrnction of each loop 
diagram, the renormalization of the composite field can be performed in a simple manner, 
as we are going to discnss later on. 

Let ns focns now on the zero modes in Nenmann bonndary conditions whenever ji = 
0(^2 7 ^ 0), = 0(ji 7 ^ 0) and ji = ^2 = 0 present in the tensor First, note that the 

zero mode composite basis fnnctions are the same, namely, = cos(ajz) and the prodnct 
involving a zero mode becomes UQ^j(z) = = ^cos(ajz). From the definition Sqj^j = 

ul)^\z)u^l\z)u^^\z)dz we read off the valnes = -^[SU + ^ 2 ) + S{j - ja)]- 

Analogonsly, it is not difficnlt to prove that = 6 {j). 

We can employ the framework jnst developed to compnte arbitrary loop diagrams in- 
clnding mnltiplicatively renormalizable composite vertex operators. At the present moment 
we shall restrict onrselves in getting the expression corresponding to the one-loop graph for 
It is represented by 



(A^ + 2) 
6 


E 

ii j 2 >o 


c{t) A(t) 




(39) 


Keep in mind that we wish to express this diagram into a nnified form. Therefore, we choose 
A A 2 7 ^ 0. Let ns start explicitly with Nenmann bonndary conditions. The snmmations split 
in the form: 


E 5.'AuAj.'L / d'‘-'qG„{,q + k,]i)G^{<,,h)^ / d"-"? LSooSSCo)? + fc, 0)G„(,, 0) 


ii j 2 >o 


OO OO 

-E'S.‘A.o4«G'o(9 + 0) + E^aojASAo)? + k, 0)G„(q,j,) + E 4hu. 

ii=i i2=i iij2=i 


^‘^jij2j^o(<I + k,ji)Go{q,j2) 


(40) 
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Using explicitly the tensor components we first obtain {i± = U ± * 2 ): 

00 n ~ 00 

ii,j2>o h=i 

00 

[h(i+ + ji) + h(i+ — ji) + 6{i- + ji) + 6{i- — ji)]Go (9 + k,ji)Go{q, 0 ) + 2'^^[6{j — ^ 2 ) 

72=1 

+ ^ 2 )] [h(i+ + ^ 2 ) + (5(*+ — ^ 2 ) + h(i_ + j 2 ) + S{i- — j2)]Go{q + k, 0 )Go(g, ^ 2 ) 

00 

+ ^ [<^(4 + ji + h) + 5(4 - ji + 4) + 5(4 + ji - J 2 ) + 5(4 - ji - 4) + 5(f_ + 4 

71,72=1 

+4) + 5(i_ - ji + 4) + 5(f_ + 4 - 4 ) + 5(7_ - 4 - 4 )] [5(i + 4 + 4 ) + 5(j - 4 + 4 ) 
+5(i + 4 - 4 ) + 5(j -4 -4)]Go(g +A:,4 )Go( 9,4) +45(4[5(4) + 5(i_)]G'o(g,0) 

X (jo(9 + k, 0)] . (41) 

In the first term of the first (second) single summation we perform the change 4 —4 (4 ^ 

—4)- This extends the range of the summation to the interval (—cxd, cxd) but precluding the 
zero mode. We can include it in order to write the hrst summation as: 

OD 

X][5(i - 4) + 5(i + 4)] [5(4 + 4) + 5(4 - 4) + 5(f_ + 4 ) + 5(i_ - ji)]Go{q + k,ji) 

jl=l 

00 

Go{q, 0) = S{j + 4) [5(4 + ii) + 5(4 ~ Ji) + 5(f- + ji) + 5(i_ — ji)]Go{q + k,ji) 

7l=-oo 

xGo{q, 0) — 2h(4[5(4) + ^iG)]Go{q + k, 0)Go{q, 0), (42) 

and similarly for the second single summation. The double summation can be manipulated 
using analogous steps and leads to 

OO 

[5(4 + ji + 4) + 5(4 “ ii + 4) + 5(4 + 4 “ 4) + 5(4 “ 4 “ 4) + 5(f- + 4 + 4) 

71,72=1 

+5(f- - 4 + 4) + 5 ( 7 - + 4 - 4) + 5(i_ - 4 - 4 )] [5(j + 4 + 4 ) + 5(i - 4 + 4 ) + 5(j 

CXD 

+ji — 4) + 5(j — 4 “ 4)]f^o(? + k,ji)Go{q,j2) = ^ [5(4 + ii + 4) + 5(i+ — 4 + 4) 

7l,72=-oo 

+5(4 + 4 - 4) + 5(4 - 4 - 4) + 5(f_ + 4 + 4 ) + 5 ( 7 _ - 4 + 4 ) + 5 (i_ + 4 - 4 ) 

CxD 

+^(^- “ ii “ J 2 )]f^(j + ji + j 2 )Go(g + A:, ji)Go (97 J 2 ) — 2 ^ + ji) + 

jl=-oo 

OO 

~ji) + 5(f- + ji) + 5(f_ — ii)]Go(? + k, ji)Go{q, 0 ) — 2 (5(j + 4) [5(4 + 4) + 5(4 

72=-oo 

~4) + 5(f- + 4) + 5(f- — j2)]Go{q + A:, 0)Go(g, 4) + 45(4 [5(4) + 5(f-)]G'o(q' + k, 0) 
xGo(?,0)- (43) 
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Replacing the value of each expression inside the combination of the diagram, the single 
summations and the independent terms cancel among each other and we obtain: 

— Q —•^2)[<^(*+ + ji + J2) + <^(*+ — ji + 32) + + ii 

ii j 2 =-oo 

—32) + — 3 i ~ 32) + + ji + 32) + — ji + ^2) + + ji — 32) + — ji 

—32)] j d‘^ ^<lGo{q + k,ji)Go{q,j2)- ( 44 ) 

Recall that we wish to express the one-loop diagram involving at maximum a single sum¬ 
mation. Performing the summations explictly, we hnd 



(A^ + 2)l 

6 2 L 


+ j) + 6{i+ - j) + + j) + - j))l2{k,j,a, hq) 


, J (i, G- 3 G+ 3 ~ \ . J (1 G+ 3 G - 3 ~ \ . 3 (1 G- 3 i-+j ~ . 
+22 — 5 cr, /io) + Gyk^ —-—, —-—, a , po) + hyk, —-—, —-—, a , /xq) 


+l 2 [k, ——, —-—, a, po) 


(45) 


2 2 

where the integrals l 2 {{k,i, a, ^q) and l 2 {k,i, j, a, are dehned in Eqs. 

For DBC the situation is even simpler since we do not have to deal with zero modes. We 
can write 


-- — Q —^ X] ^qGo{q + k,ji)Go{q,j2)- (46) 

iij 2 =i 

Using the tensor components already derived, we can write the summation as 

00 p CO 

y-ii2iii2y'i]2j J ‘G ^?Go (9 + ^,ii)G'o(g,i2) = [d{G + 3 i + 32) + GG — 3 i 

31 , 32=1 31 , 32=1 

- 32 ) - 5(f- + 3i + 32 ) ji - is) -d{i+- ji + j 2 ) -5{i++ ji - 32 ) + - ji 

+ 32 ) + d{i- + ji - 32W{3 - h + 32 ) + G3 + 3i - 32 ) - d{j - ji - j 2 ) - 6{j + ji + J 2 )] 

X J d'^ ^QGo{q + k,ji)Go{q,j2)- (47) 



Now, in the second bracket perform the change j 2 —j 2 in the hrst and fourth terms. This 
has the effect to produce a global factor of (—1) in each coefficient of those two terms. The 
net result can be written in the form 


CX> 

E o(l) o(l) 
‘^hi23l32‘^3l323 

31,32=1 


d'^ ^qGoiq + kJi)Go{q,j2) 


CO CO 

GJ + 3i + 32)[GG + 3i 

ii=ii2=-ooj2^o 


+ 32 ) + GG- 3i - 32 ) -d{i-+ ji + is) - GG - ji - is) - 5(4 - 3i + 32 ) - 5(4 + ii - is) 
+5(ii--ii+is) + 5(x_-Fil-is)] j d'^"^gG'o(gfc,ii)Go(g,is). (48) 
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By performing similar transformations in some terms involving ji{^ —ji), we can extend 
the summation over ji for negative values as well. The simplicity here is that no zero mode 
survives when we include the values ji,j 2 = 0 in the summations. The graph then reads: 

^^6 + ji + h)[S{i++ ji + j2) - S{i+- ji + j2) - S{i++ ji 

ii j2=-oo 

-j2) + S{i+- ji - j2) -S{i-+ ji + j2) + 6{i- - ji + j2) + S{i-+ ji - J 2 ) - S{i-- ji 
-J2)] j (f~^qGo{q + k,ii)Go{q,j 2 ). (49) 

Comparing this with Eq. fHT)) . DBG (r = —1) and NBG (r = 1) can be unihed very easily. 
When this diagram is written in terms of the integrals J 2 (A:, z, d,/ iq) and J 2 (fc, z, d,/ iq), the 
unihed result is 



{N + 2 )t 
6 2 L 


+ j) + <5(4 - j) + r(5(z_ + j) + (5(z_ - j)))l 2 {k, j, a, po) 


4-i 4 + i ~ N T,i G~j ~ M f/, ^- + J ~ N 

+^[h{k, ——, ^ , a, /io) + h{k, ——, —-—, a, po)] + 4(fc, —^—, cx, po) 


+4(^, ——, cr, /zo) 


( 50 ) 


This can be rewritten in a more elegant, compact notation as: 


{N + 2 )t ^ 


6 2 


-^ 02 " 5(4 - ^j)/2(A:,4,cr,;[Zo) + 


4 + 13 j ia - f3j 


a, /io) 


a,B=-l 


( 61 ) 


This framework has the virtue of expressing all integrals in terms of inhnite sums which 
makes it simple the comparison with the results from periodic and antiperiodic boundary 
conditions. Needless to say, all the massless integrals follow from the substitution /zq = 0 in 
the integrals above and in the higher-loop contributions, which will be analyzed next. 


3. Higher loop diagrams for and 

Our goal now is just to write down the various graphs in terms of integrals which resemble 
those from periodic and antiperiodic boundary conditions, possibly with additional “non¬ 
diagonal” terms. We are going to consider some two-loop graphs, although just a smaller 
subset of them will be necessary to our computation of the critical exponents within the 


21 






















present technique. We will analyze only one three-loop graph for the two-point function 
which is the only one needed for our purposes. We will save the construction of other non¬ 
trivial diagrams (external momentum dependent) including mass insertions for later when 
we will discuss the renormalization of this theory. 

Start with the two point function for instance. We can classify the diagrams in trivial 
and nontrivial contributions. The trivial contributions will be generically called “tadpole 
diagrams” and are going to be discussed hrst. These contributions are displayed in the 
two-loop case for the sake of comparison with the one-loop case but we will not need them 
in the present computation. 

Of course, we can start from scratch with two four-index hnite-size tensors and perform 
the apropriate contractions between them in pretty much the same way we did in the one- 
loop case. The result including both boundary conditions is: 

~ % % 

(5(z_)D2(d, Ho)Di{a, /xq) + r/2(0, y, y, d, Ho)Di{a, Ho) -F r/2(0, z_, 

d,;Uo)I)i(y,d,/ro) + rh(4)T)2(d,/io)T»i(d,/io) + / 2 ( 0 , y, y, d,/io)T'i(d,/ xq) + h{0,i+, 

OD . . . . 

^ z ^ r ~ z z ~ ^ I 

5‘,^o)-Di(^,a,/xo) + 22 ['^ 2 ( 0 , ^ ^— ^5‘,/zo)-Di(/,5‘,/io) + r/2(0, ^ + Z, ^ a, 

l=—oo 



a, fio) , 


(52) 


where 


= (53) 

L = — 00 

The nontrivial two-loop graph of the two-point function can be obtained most easily 
by taking the four-point graph, making the identihcation ii = is = j, introducing a new 
propagator with momentum q 2 , mode j, performing a summation over j and an integrating 
over d'^~^q 2 - We hnd: 





(d(i-) -F Td(i+))l3(k, ii, a, /io) + ^hik, 


!± 

2 ’ 


i- 

T’ 


(T,/io) + Srls^k 



4 
2 ’ 


( 54 ) 
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where, 




^qid!^ ^q2 


J2=-0O 


[qI + 5-2 j| + /ig] [(gi + q2 + kf + 5-2 (ji + ja + jf + ho] 


X 


fe? + + hg] ’ 

CO 

'2 


(55a) 


Ig^ 




X 


i(9i+®+*p+*^((+jp+/ia' 

The nontrivial three-loop graph contributing to the two-point function can be determined 
similarly. The complete result is: 




{N + 2){N + 8) 
27 


+ T 6 {i+))l 5 {k, ii, a, po) + h{k, y, y, a, po) + rhik, y, 
y, d, Po) + y, y, * 2 , d, po) + ^^{k, y, y, Zs, d, Po) + 2/5 (/c, y, y, Zi, d, po) + 2 r /5 

y) y) *1) d, po) + h{k,ii, y, y,d,po) + y, y,d,po) + Z5(h,zi, y, yjd, po) 


- 1 -t 75 (A:, Zi, y, d, Po) + 75 (h, y, y, Zi, d, Po) y, y, —Z 2 , d, Po) 


(56) 


where the integrals above are dehned by (the arguments (fc,/, m, rz, a, po) were suppressed 
in the integrals with three external quasi-momenta indices for sake of simplicity) 

d‘^~ ^ qi d'^~ ^ g 2 ^ gs 


n 

h{k,i,a,iio) = Y j 7^ 


J1J2 J3=-00 


+ d^ii^ + ^2] [g2 ^ ~2j2 ^ ^2] [g2 ^ ~ 2g2 ^ ^2] 


X 


[(gi + g2 + ky + (T 2 (ji -K j2 + z)2 -h ho] [(91 + 93 + + Cr^(jl + js + *)^ + ho] ’ 

d^~^ qid^~^ q2d^~^ q'i 


(57a) 


n 

h{k,i,j,d,qLo) = Y / 


Jlj2 = -00 


[g^^ + d2z2 + //2] [g2 ^ ^2g2 ^ ^2] [g2 ^ ~2j2 ^ ^2] 


X 


[(gi + 92 + ky + a^{ji + jy + pg] [(gi + gs + ky + a‘^{j2 + jy + ho] ’ 
d‘^~^ql2{q + k, l,j + m)/2(g + k,i + n) 


n 

& = yE/ 


^ ~ ^ f d^-^qhiq + k, g 2,92 + ji + ^^ 2(9 + k, 92 + m, 92 -F ji + n) 

jl J 2 = -00 ■ 

°° r xd— 1 , 

4 = 5 E 

jl J 2 =-oo ■ 

°° /■ ^d— 1 , 

4 = 5 E 

9 l J 2 = -00 ' 


g 2 + ^ 29-2 + ;x 2 

d^-^qhjq + k, 92,92 + ji + 072 (g -F fc, 91 -F 92 + m, 92 -F n) 

92 + (729-2 + 9^2 

^*^"^9^2( 9 + 7 , ji + 7 rn)i2{q + /c,92, 292 -F ji + n) 
g 2 + ( 729-2 + p 2 


(57b) 

(57c) 

(57d) 

(57e) 

(57f) 
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Let us turn our attention to the four-point vertex part. The simplest contribution at 
two-loops is the double bubble. Using the fusion of two one-loop diagrams of the four-point 
function results in the following unihed expression {l 2 {k, d, po) = hik, z), l 2 {k, i,j, a, po) = 

{N^ + QN + 20). 


XXX^ 
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a\ [5(z+ -+- z+) -F r(5(z_ + z+) 6{i+ - z+) -F rh(z_ - z +)]/2 (fc, i+) 


+ [r5(z+ + z_) -f 5(z_ + z_) -F rh(i+ - z_) 5{i_ - iJ)\ll{k, i_) 2r(/2(A:, i+) + hik, i+)) 

X h(k, + 2T(h(k,'i.) + h(k,t.))h(k, ^=4^) + 2{h(k:i.) 

+ h(k, U))Uk. + 2(h(k, ?+) + h(k, i.))h(k, 

oo 

+ 2 ^^ (l2{k,j,j + i+) h{k,j + 


+ rl2{k,J + 


j=-oo 


z ,, . . . sf?/, ■ i-— i- . i- + i-, ~ — i+ . i_ + i+. 

+h{k,],] +i-)yi2{k,] + ^^ -) +rl2ik,j + ^>2 +-^-, 

The important point which simplihes our task is to realize that the last terms like, for 

OO _ _ . - . - 

instance, Y1 h{k,j,j + i+)l 2 {k,j + , j ) are 0 (e°) and do not contribute to 

j=-oo 

the singular part of this diagram. These terms can be ignored. 

Consider the evaluation of the nontrivial two-loop contribution of the four-point vertex 
function. Following the same line of thought we obtain {P = ki + ^ 2 ) 

(5A^ + 22) 


(58) 




27 


a< [5(4 - 4 ) -h t 5(4 - z_)] 14(4 4 , 4 , Zg) -F [5(U - z_) r5(z_ - 4 )] 


x/4(P, 4 ,z_,Z3 ) -f [5(4 -F4) -Fr5(4 + i-)]h{P,k 3 ,i+,-i 3 ) + [5(z_ + i-) + T6{i- +4)] x 

T/n 7 ■ -N f/n 7 4 + *- 4 “ f/n 7 i-+ i- i-— i- i-— i+s 

4(4 4, -^3) + 4(4 k3, —-—, —-—, —^) + r4(4 ^ 3 , —^^—, —^) 

T/T^i *+“*+*+ + *+ 4 “ , f / 7-, 7 i- — i+i- + i+i- — i-., f/n7 4 

+''■ 4(4 4, —^^^—) + 4(4 4, —^^^—) + ''■ 4(4 4, Y 

4 . ^ 3 . 4 ^. + 4(P. * 3 . ^. 4^4 

+y) + 7/4 (P,*s. 4^' 4^’ 4^* + 2p4(n *:s.7+. 4^' 4^' 

. z_-hz- z--4n f 7 ■ 4~4 4 + *-\ f/n 7 • z_-4 + 

*-7 —^—7 —+ ^4(4 47 *+7 —^—7 —^—) + 4(4 4, *-7 —^—7 —^—) + 

h, i+, 4 ^, 4 ^) + pp. k,, i_, 4^, 4^) + '■•3,7+, 4^. y 

-y) + /4(P,fc3,7-,4^>4^']}’ 


where the objects h{k,k\i,i){= h{k,k\i, j]a, ^o)),h{k,k' h{k, k',i, j,l; a, fio)) 
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and l 4 ^{k, k' ,i, j,l){= l 4 {k,k',i,j,l;a, fio)) are defined, respectively, by 


h{k,k',i,j) = a'^ 


l,m=—oo '' 
1 


X 


[{qi -q 2 + k'Y + cr2(/ - m - j )2 + ^l] 


[ql + - kf + a\l - if + + /^o] 

(60a) 


21 ’ 


oo « 

h{k,k\i,jf)=a‘^^ J 




m=—oo 
1 


[qf + aH'^ + /ig] [(gi - kf + cr2j2 + /ig] [g 2 ^ 2^2 ^ 2 ] 


[(gi - g2 + k' f + cr2(/ - m)2 + p2] ’ 

CX) „ 

h{k,k',i,j,l) =d^ Y1 

m 


(60b) 


^gid'^ ^q2 


m=—oo 
1 


X 


2 + (T2m2 + Po][(gi “ kf + a‘^{m- if + qf][q2 + + if] 

(60c) 


[(gi - g 2 + k' f + cr2(/ - m)2 + p2] • 

Note that the pole in e coming from the computation of the summation is absent in the 
integral / 4 . Since we are only interested in the singular part of the above diagram, we can 
simply neglect the contribution coming from this integral. 

We are now left with the two-loop diagrams of the vertex part 
r‘^^’^^(/ci, ^ 2 ; <5, * 2 , j,/^o, d). Sticking to the method above described using a similar 

simplifying notation for the integrals by omitting the arguments (/io,d), whereas perform¬ 
ing the identifications {ki) = {ki, k 2 ', Q), k = ki + k 2 , it is not difficult to find the following 
expressions {i± = f in 12 )'. 


S: 


ih) = 


r {N + 2) 
2 108 


2 


{5{i++ i) + 5{i+- j) + T5{i_+ i) + T5{i_+j))ll{k,j) 


+2rh(k, j) + h(k, 4)] + +2h(k, j) + h(k, *_)l 

+2f; (h(k,td±±2L)Uk,I i + 4) + Uk, 


l=—oo 

X i2{k^l^l i— 


(61a) 


(ki) = ^ + j) + S{i+ - j))h{k,Q,i+,i2)) + T{6{i_ + j) + S{i_ - j)) 

xlfk, Q, ^ 2 ) + Ifk, Q, + riff Q, ^^) + ifk, Q, 

i- j -i- j + i+. Y j -i+ j + i+ j -i-. . „ f ^ ■ J - 4 i “ ^ 

—)—^^—) + Th{k,Q, —-—,—-—,—-—) + 2Tl4{k,Q,j, —-—,—-—) 


I of ■ j d- i- J + 4 \ I o f (1 rt ■ + 4 J + ^ ^ ■ 3 ~ 3 ~ f 

+2l4{k,Q,j ,—-—,—-—) + 2Tl4{k,Q,j ,—-—,—-—) + 2l4{k,Q,j ,—-—,—— 
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We have now a set of integrals which are very similar to those coming from periodic 
boundary conditions (see Ref. |22j). From our discussion, getting the expression of the 
diagrams is the laborious part of the method. Since we have already discussed the periodic 
integrals, only minor modihcations are necessary to perform the computations aimed as we 
shall see next. 


III. MASSIVE FIELDS FOR DBC AND NBC IN THE EXPONENTIAL REPRE¬ 
SENTATION 


The vertex parts depend on the boundary condition, but we will suppress that dependence 
in what follows. By now it is rather clear that our unihed description in the last section leaves 
no ambiguity to dehne the renormalization algorithm simultaneously for both boundary 
conditions. 

Before discussing explicitly the renormalization algorithm, it is interesting to restrict 
ourselves to the minimal number of diagrams to perform our computation of the critical 
exponents. The argument for the massive theory regarding inhnite systems appeared recently 


25| . We will summarize the steps here and use the results directly. Start by dehning the 


three-loop bare mass /tq in terms of the the tree-level bare mass as = 0, j, a, /xo, Aq) = 

AS = = o.i , C7, /iO) -^o) “ Then, inverting this equation to get the tree-level //q 

in terms of the three-loop bare mass /Iq and all the diagrams computed at zero external 
momenta has the virtue of eliminating all tadpole diagrams for they do not depend on 
the external momentum. Next, expanding /ig = /no(/ho,Ao) in each primitively divergent 
vertex part at hrst order in Aq eliminates all the remaining graphs containing one-loop mass 
insertions. Here is one of the most nontrivial results and we are going to discuss this topic 
carefully in what follows. 

Recall that the two-point function always includes a factor Si^i^ = A(zi — * 2 ) in arbitrary 
loop diagrams. The set of steps just described permits one to write the tree-level bare mass 
in terms of the three-loop bare mass up to 0(Ao) through 




Ao 




( 62 ) 


Notice that when we perform this expansion in any integral, a complication occurs: the 
tadpole graph has hnite size (kind of “internal symmetry”) indices, say l,m, which will be 
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attached to a propagator of a certain internal line of an arbitrary graph. So these indices 
should be contracted with the set of indices taking place naturally in the diagram. Consider 
the two-loop diagram of the four-point function with a “mass insertion”: 


r) ^ {N + 8){N + 2) yy 

il,j2 ,i3 j'4>0 

Go{Qi + k, j3)Go{qi, j2)Go{q2, ji)- 


(f ^qi(f ^g2Go(gi + fc, ji) 

(63) 


(The condensed nomenclature = ^[iii 2 jij‘ 2 \T is sometimes useful.) In this expression 

we can set /xq = /fo since the correction will produce three-loop terms which are not going 
to be useful to our analysis of the four-point function up to two-loop order and is going to 
be neglected. 

Let us take a look in Eq. (IT^ . There the propagator is dehned as Go{q,j) = Go{q,j, ^q). 
When we replace the mass po by fio into the propagator, the left of the graph carries 
index j from the original propagator before the insertion. Due to the internal character of 
the propagator which carries hnite-size symmetry index, the coupling of the left hand side 
with the inserted mass tadpole is implemented as follows: the index j must be free in the 
resulting expression and all the the other indices appearing there should be contracted with 
summations. This is equivalent to the following simple rule 


Go{q,j,^io) = Go{q,j,fio) 1 - X] 


An 


l>0 


2[g2 + 


XI 


n -1 




Go{<l,j^f^o) + ^ G'o(g,i, po)G'o(g,/,// q) _Q_ 


l>0 




(64) 


Now, replace this result in Eq. ([T9]). The one-loop diagram has two propagators which 
contribute the same amount to the correction due to the mass insertion (after some redeh- 
nitions). Upon substitution of Eq. (|8|) for the tadpole, we hnd 


o 


(iV + 8) 
9 


E / d‘^-'qGo(q + k,h,Mo)Go(q,h,/Io) 




(N + 8)(N + 2) ~(r) 
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c(r) ^(t) 

j2j3i3U 


(f ^qi(f ^g2G'o(gi + ^,ii,ho) 


Jl J'2J3 ,j4>0 

Go{qi -f A:, j3, fio)Go{qi, j2, /fo)G'o(g2, j 4 i ho)- 


( 66 ) 


Then, we conclude that the two-loop correction term due to the expansion of the internal 
propagators in the one-loop diagram of the four-point function produces the mass insertion 
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whose value is equal exactly to that from Eq. 
important graphical identity: 


This statement is equivalent to the 


)0( 1/^0 - Uo+^ )Q( 


I Mo • 


( 66 ) 


By following the same prescription, all the primitively divergent vertex parts with mass 
(tadpole) insertion are eliminated. For instance, the reader can check from what we have 
been discussing that the following diagrammatic identity also holds 

3An 




I Mo 


(67) 


It is easy to conclude that similar diagrammatic identities are valid for the composite op¬ 
erator due to its similarity with Despite the apparently complicated tensorial 

structure of this hnite-size held-theoretic formulation, the mass insertions are cancelled for 
arbitrary external quasi-momenta modes as shown above in a simple manner. This cancel¬ 
lation takes place since different orders in Aq in the perturbative expression have different 
signs, and the combinatorial factors coming from the mass insertion corrections in each di¬ 
agram exactly matches those from the perturbation expansion of the original diagram (one 
power of Aq higher) in fiQ. This suffices to prove the claim explicitly at the perturbative 
order considered in the present paper. See 2^ for more details. 

The resulting vertex parts now depend exclusively on /io, Aq. Their diagrammatic expan¬ 
sion in terms of the reduced number of diagrams are: 

A^ 


= k'^ + P-l- ~ 

l/io 


\k=0,p.Q 


+ 




Ah 


\k=0,iio 


^hLisiAki^ cr, po, Ao) — XqS, 


\2 

-^0 




4 V 


)CX],h‘ + ^{([-^ 


lAiO 


(68a) 


(ki) -|- bperms. j + 2([>OCX 


rS,%(fci, k 2 ; Qs, a, /io, Ao) = rTlf-^ + T 


{ki + k 2 ) + 2perms. 

Mo 

(2,1) y.(2,l) 

n,*2,1=0 *i,*2,i=*i-l-*2’ 


MO 

(68b) 


^hliljiki, k 2 ] Qs, cr, /io, Ao) — 


3A, 


+6An 


{ki,k2',Q3). 

J Mo 


_ (^1 + k2] Qs) + 3Aq 




(68c) 

{ki -|- k2] Qs) 

(68d) 

for the vertex T^^^ since after using 


We choose not to display explicitly the tensor 
the value of this tensor, the two-point function depends upon only one external quasi¬ 
momentum (mode). Note that Eq. f|68cp takes into account that each diagram appears in 























the combination mentioned before in Eq. (1381) . We appplied a similar simplifying notation 
for the vertex part and choose not to show the dependence on the tensor in the 

zeroth order diagram because of that implicit combination. (We could have done the same 
to the vertex part provided we include the zero order graph in the diagrammatic 

expansion.) We omit hereafter the lower indices (“internal” hnite-size modes) whenever 
referring to an arbitrary vertex part but keep then on the argument of that vertex. 

The vertex part j, a, fiQ, Xq) = T^‘^\k,j,a,jlo,Xo) — has the advantage of 

not possessing the tree-level term and looks like the two-point vertex part from the 
bulk case. It has a logarithmic divergence as well as ji, a, flo, Xq). In our specihc 

case of a given order in the coupling constant expansion, we want to go up to three- 
loop order in the expansion of Dehne the renormalized quantity f^^(/c, j, d, po, Aq) = 
, cr, po! Ao). The introduction of the normalization function produces the can¬ 
cellation of the logarithmic divergence aforementioned. After that, we dehne a renormalized 
hnite mass at three-loop order as On the other hand, the four-point vertex part 

starts proportional to the bare coupling constant, so we can dehne a renormalized coupling 
constant up to two-loop level by writing g = Z'^Xq. First, write = l + gzi + g‘^Z 2 + g^z^ +... 
(where Zn are divergent quantities). From the absence of linear terms in the coupling con- 
stant perturbative expansion of E}^ , we immediately obtain zi = 0 . At third order, we can 
express it entirely in terms of renormalized quantities, f^^(/c, j, d, p, g) = f^^(A;, j, d, po, Aq). 
The same happens to the other vertex part, since up to two-loop order the renormalized 
object dehned by {ki, ji,a, g.., g) = Z^E(‘^^(/c, jj, d, po, Aq) is automatically given in terms 
only of renormalized mass and coupling constant, as the reader can easily verify. The same 
argument carries through to the renormalized composed vertex part obtained from the log¬ 
arithmic divergent bare one via E^’^^(fc, j, d, p, (?) = Z 02 Z 0 E‘^^)(A:, d, po, Aq). It requires 
another normalization function Z ^2 and no longer depends on bare quantities at this order. 

In general, we can go to arbitrary loop order and get rid of all the bare quantities in 
the dehnition of generic vertex parts which are renormalized multiplicatively. Multiplicative 
renormalizability amounts to say that an arbitrary vertex part E('^’^)((L, M) ^ ( 0 , 2 )) in¬ 
cluding composite operators can be renormalized through the functions Z ^^'^, Z ^^2 such that 
the vertex parts dehned by 


p(L,M) 

^ R 


{Pni ^ni 


Qn' I jn': 9 1 h) 


(A A( AVrli'">(p„, i„, Qn'Jn; A„, Ao, A), 


(69) 
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are automatically finite. (In the argument above for {L,M) = (2,0).) 

Here, g, g and A are the renormalized coupling constant, mass and cutoff, respectively. 
The argument works well for all kinds of regulators, which are used to express inhnite 
(divergent) quantities in terms of functions of these parameters (regulators). We employ the 
cutoff when deriving differential equations satished by renormalized vertex parts. From the 
explicit operational viewpoint, however, we shall focus here in dimensional regularization, so 
that divergences are written in terms of poles in e = 4—d (oc e~\ where I is an integer positive 
number, usually the number of loops), where d is the space dimension of the system. The 
symbol n labels the external momenta pn as well as the mode of the external quasi-momenta 
in associated to the external legs (n = 1,...,L) of the helds. The label n' is connected 
to the external momenta Qn' and the mode jn> of the quasi-momentum corresponding to 
the external legs {n' = 1,..., M) of the composite helds {(jP‘ insertions) in an arbitrary IP I 
diagram. 

The normalization conditions on the primitively divergent vertex parts for the massive 
theory are chosen at zero external momenta and nonzero quasi-momenta. From the structure 
of the diagrams, it is obvious that the condition on the external momenta are not sufficient 
to simplify our task. In fact, there are other renormalization schemes, for instance minimal 
subtraction, which do not require any hxation of the external momenta. In order to give 
a simpler prescription useful for all other renormalization methods required from an ab 
initio formulation, we choose the external quasi-momenta to be set at especial values. The 
symmetry point here is dehned at external quasi-momenta from the external legs chosen to 
be equal and set to the value i (arbitrary nonvanishing positive integer). The hniteness of 
the theory can be achieved through the conditions 


= 0,i,g,p) = 


dT^R\k,i,g,p,) 


dk^ 


k^=0 


= F 


r);\^z = o,p = *,^,p)^Fg) 


R ^ 
d2,l). 


SP 


= Sdfif, 


Tr’ ik = 0,ii=i2 = i,Q = 0,j, g,p) = T 


— t(2,1) 

R 


SP 


3r 

Y' 


(70a) 

(70b) 

(70c) 

(70d) 


We can now discuss the renormalization group invariance of arbitrary renormalized vertex 
functions. At the critical dimension d = 4 utilize the cutoff as the regulator. The inhnite 
cutoff limit in the integrals appearing in each vertex part multiplicatively renormalized is 
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taken solely after the renormalization prescription is established. Apply the derivative 

0^0 

over the bare vertex part ((A^, M) ^ (0, 2)) at hxed Aq, A in order to obtain the vertex 

function Qn'0; Aq,/ io, A) at zero inserted momentum. Next, rewrite the 

remaining bare vertex parts in terms of the renormalized ones. After similar manipulations 
performed in the case of periodic and antiperiodic boundary conditions in 22 | but utilizing 
the above normalization conditions, one learns that the renormalization group invariance of 
the renormalized vertex parts in different renormalized mass scales is substantiated in the 
following form of the Callan-Symanzik equation: 


(71) 


dlnZ 




where = p.|^, qJ’ and 7 ^ 2 ^ = —jj —Even for z 7 ^ 0 in the present 

unihed framework, this form resembles very much the situation for PBC. 

When the held is not at the critical dimension, the annoying dimensionful aspect of 
the coupling constant can be circumvented by writing the bare (renormalized) coupling 
constant in terms of a genuine dimensionless bare (renormalized) coupling constant Uq {u) 
as Aq = (g = g^u). The how function g) can be written in terms Gell-Mann-Low 

function dehned by [f3{g,g)]GL = —^g + I3{g,g)- Using the Gell-Mann-Low function into 
the CS equation, along with the dimensionful quantities dehned in terms of dimensionless 
amounts, we hnd that [(3{g, g)]GL-§^ = turns out that resulting renormalization- 

group picture involves only the dimensionless renormalized coupling constant and possesses 
a well dehned scaling limit 


dlnZ 


(t) 
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30j. The Gallan-Symanzik equation now reads 
L 

~ ' 2 ' 




dg du 

/^^(2 - 7j,^^)rR’^^^^(Pn,G,Qn',C;0>U,p) , 


(72) 








7^2^ (u) 


dlnu, 


7 )' 


= —e 


du 


= /?(")(«) ' 


= /3(")(n) 



(73a) 

(73b) 

(73c) 
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where 











Another function that will be useful to our purposes and utilizes the definition = 


is written as 




= /?( 


t) 


d lnZ^J ' 

du 


(74) 


The solution of the Callan-Symanzik equation has already been described in the ultraviolet 
scaling regime. The right-hand side becomes negligible at the scaling limit, ^ cxo j^. 28|. 
showing that the solution admits a ultraviolet hxed point. We will just make use of these 
results here and prove that the dehnitions above are sufficient to our computation of the 
hxed point and its aftermath in the evaluation of critical indices. 

Let us start the discussion of the asymptotic limits for both boundary conditions. The 
basic objects to be computed are the one-loop integrals hik, z, d, /r) and i,j, a, which 
belong to the one-loop graph of the four-point function. The resemblence of these integrals 
for DBC, NBC, ABC and BBC will per mit a more economical approach to this topic and 
the reader is advised to consult Ref. (22| to grasp many details. We prefer to omit them 
herein. 

The integral l 2 {k,i, a, fi) is identical in the form to l'^\k,i,a, for periodic boundary 


condition (r = 0 in Eq. (10) from 22j). First, we factor out the mass fi from the integral, 
rescale all momenta and dehne h = A Second, we introduce a Feynman parameter and solve 
the integral over the momenta. Third, we perform the summation using the representation 
from Ref. 2^. Fourth, transform the argument of the resulting (factorial) T function from 
{d — 1) to d. Then, divide the integral by the area of the unit sphere at d dimensions. These 
sets of steps lead to the following result (see also Appendix A) 


h{k,i;r) = 


l2{k,i;a,fi) fi' 


Si 


-—/ dx[x{l — x){k'^ + + 1] 2 


r(2-0q(fc,i;f)), 


where 


Fa{k,r,r)=r y dxfi^^\^xi,h{k,i,r) 

00 

fa{a,b) = 4'y~^cog(27rma)(—p^)" K^_^{2nmh), 

m=l 

h{k,i,f) = f~^\Jx{l — x){k? -|- -\- 1, 


(75) 

(76a) 

(76b) 

(76c) 


where K,,{x) is the modihed Bessel function of the second kind. 
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Using a similar chain of reasoning, it is not difficult to find that the other integral required 
reads: 

/2(fc,i, j;f) = ^ dx{r‘^[xi‘^ + (1 ~ — a;) + l) (77) 


Recall that in order to compute the one-loop diagram, owing to our use of nor¬ 
malization conditions in this massive framework, we need the previous integrals eval¬ 
uated at vanishing external momenta {k = 0). Not only this: specihcally, we need 
12 ( 0 , 2f, r), 12 ( 0 ,0, r), 12 ( 0 ,0, 0; f), 12 ( 0 , h 7) and / 2 (/i:, 2z, 0; r). In that case, setting A; = 0 
and e = 0 in the 0(e°) terms from l 2 {k, 2z; r) we hnd: 

/2(0,2i,r) = / ln[4:x{l - + l]dx + ^Fo{0,i,r)\ (78) 

It is also simple to demonstrate that J2(0,0,0;r) = / 2 (A:, 2z, 0; r) = and 

l 2 {k,i,i;f) = 2 ^ 1 +fii^ ■ Replacing these results in the expression of the one-loop four-point 
diagram for the particular choices of external momenta and quasi-momenta, we can write 






9e 


1 - f + 7) 


CR)(z; 7) = iFo(0, 2q f) + lFo(0,0, 7) -1- f + 


2rT 


+ 


2f 




arcsinh{ri) — 1 


(79a) 


(79b) 


In the above equation r oc (j)”^ where ^(~ /i~^) is the bulk (inhnite system) correlation 
length. We will analyze this correction term in order to study its asymptotical behavior for 
large as well as small values of L. 

Let us start by the hrst scaling regime, namely j —)■ cxd (7 ^ 0) limit. It is easy to 
demonstrate that i) the linear terms in r go to zero, ii) Fo(0,0,7 —)■ 0) ^ 0 [22| and in) 
the last term also vanishes. It is not difficult to demonstrate that Fo(0,2h7 0) tends 

to zero as follows. From Eq. fl76cp . limh{0,2i,r) —)■ r~^, and we can solve the integral by 

f^O 

writing explicitly the integrand in terms of a summation involving Bessel functions. The 
coefficient of each term in the summation is equal to sm(47rmz) which is zero for integer i, 
therefore proving the assertion. Then, the hnite-size correction reduces to the bulk result 
(^(t)(^. ^ ^ ^ whenever | cx). 

Next we consider the limit 0 (7 —)■ cxd). The several terms inside the correction func¬ 
tion are going to be examined separately. Note that lim^^ — lim ^ —)■ It, 

^ ^ n J r^oo3[l+Ul+4^] 32’ 32’ 
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which are convergent and r-independent. In addition, lim — ^ arcsinhiri) — 1 

— |[/n( 2 fz) — 1 ], and we have a logarithmic divergence in the desired limit. 

We now focus on the contributions from FQ{0,2i,r) and id)(0,0,f) in this limit. The 
function 

CXD 

fi{2iXy h(k — 0, 2z, f) — 4^cos(47rmia:) A'o(27rm/i(fc — 0, 2i,f), 


( 80 ) 


m=l 

can be better understood through the identity 

OO 1 ' ^ ' 

Ko{nx)cos{nxt) = - 7 +/n(- 

n=\ 

'E 


2 l| 


dvr' 


+ 


TT 


2 £^/^TI 2 


TT 

2 


n=l 


1 1 

OO 

1 1 

'^x2 + ( 2 n 7 r + tx)2 2 n 7 r 

^ 2 ^ 

n=l 

'^£2 + (2n7r — tx)2 2n7r 


( 81 ) 


where 7 = 0.57721566... is the Euler-Mascheroni constant. It is valid for positive values 
of the variable x. Performing the the identifications x = 27rr“^A/l + 4:f‘^Fx{l — x), t = 
, we get to 

yj l-\- 4 r‘^i^x{l—x) 


fi {2ix, h{k = 0 , 2i, r) = 27 + 2ln 


^1 _j_ — x) 


2r 


+ 


E 


+ E 


\/l + Ar'^Fx 
1 1 
n. 


^ ^ \/+ 4i(i — n)x + ^ \/+ 4i(i + n)x + m? 

In particular, for i = 0 


(82) 


CXD 

/i (0, h(A; = 0, 0, 7) = 27 - 2ln{2f) + 7 + 2 ^ [ 




^1^7-2 + 


71,2 n. 


(83) 


This coincides exactly with id)(0,0,r). Taking the limit (setting r = 00 into the last term) 
we find /fmFo(0, 0,7) f — 2ln{2r). Finally, we have to analyze Fo(0, 2i, 7) which is given 

r—)-CxD 

by 


Fo(/c = 0 , 2 z, r) = 27 + 2 / dxln 


X 


-\/r + Ai{i — n)x + 


n 


a/1 + APFx{l — x) 
27 

OO „i 

+ / dx 


+ dx 


^ Jo Ly/r 2 -|-4^(^-|-ri)x + n2 


0 \/l + 4 r 272 x 

1 - 
n. 


00 „i 

+ V / dx 
(84) 


One can show that /im dxln 


yj l+ 4 r^i^x(l—x) 


2 r 


—)■ /m — 1, which is finite in this limit. The 


third term can also be shown to be finite, namely, lim —)■ ( 7 ))- 
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The summations conceal an underlying subtlety for z 7^ 0. Looking at Eq. there was 
no problem whatsoever in replacing directly r = 00 when z = 0 at the summations. That 
was the way we renormalized PBC and ABC with z = 0. Indeed, setting directly (r = 2r in 
those boundary conditions) = 0, these terms contribute zero to PBC (just as here) and 
{2ln2 — 1) for ABC. In other words, they contribute a finite constant r-independent, which 
do not have a pathological behavior in the limit r —>■ cxo. Looking at PBC, if we do not 
choose to renormalize the theory at z = 0, but choose instead all external quasi-momenta 
equal to z > 0, then the form of these summations is identical to those which appear in 
DBC and NBC here. In ABC the occurrence of this feature is a little bit different since 
one has to add the hnite term above mentioned. 

Now consider z 7^ 0 and see what happens when we set = 0 in the summations. 
Using the power series expansion from (1 -f- y)~^, keeping all terms in y and performing the 
integrals over the Feynman parameter x, we obtain 


/ 00 ^ 

Urn \ Jq dx 


\ n=l 


. r~'^-\-4i{i—n)x+n'^ 


CXD 

+ T. lodx 




2J2r"C{^n+l), 

n=l 


(85) 


where C(2ir+ 1) ^ 

p=l 

It is clear that the series diverges for z > 1. This divergent number is harmless in the 
limit f ^ oo, since i) it does not depend on f, ii) it can be eliminated by suitable choices 
of external quasi-momenta in PBC, NBC and ABC {e.g., i = A). As far as the series does 
not depend explicitly on r, it can be safely considered regular in this limit in comparison 
with poles in e. In order to eliminate the divergent term for PBC, ABC, DBC and NBC 
when we choose z 7^ 0, which is clearly an artifact of the method, we can dehne a “normal 
ordered” (regularized) correction : •= | *^"'C(2zz-|- 1) that is still 

divergent in the limit L —)■ 0 but without the nonphysical divergence associated to this 
summation. This normal ordering operation ammounts to say that if one replaces directly 
r = CX3 in the inhnite sumations Eq. flSip (or alternatively into Eqs. fl82p . fl84p ) whenever 
z 7^ 0 in analogy with which is done in the case z = 0, one shoud be careful to subtract this 
term aformentioned. That is the actual quantity that should be used in the limit of small 
values of L. 

Putting together all the terms examined in this limit, we hnd that the regularized hnite- 
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size correction function reduces to 


lim : r) —(2r + 5) — ln{2r) + 

f^oo Qi 3 

The finite-size correction is clearly boundary condition dependent. From this expression 
we can obtain two distinct behaviors in this limit whenever the system is not located at the 
bulk critical temperature (f 7^ 0). For DBC{t = —1), the divergence is purely logarithimic 
like in ABC, whereas the dominant term for NBC{t = 1) diverges linearly with r just 
like in BBC. When these terms become comparable to the dimensional pole in e through 
decreasing values of L the dimensional crossover takes place. 

Let us turn now our attention to the two- and three-loop diagrams. They are required 
to getting the critical exponents r] and u perturbatively in the safe scaling region discussed 
above where L is not so small and e-expansion results are valid. We will just write down the 
solution of the integrals. The interested reader should consult the details in Appendix A. 

Consider the two-loop diagrams of the four-point vertex function. They consist of “diag¬ 
onal terms” which are responsible for the leading singularities in e as well as “nondiagonal 
terms”. The double bubble can be easily obtained from the one-loop diagram just discussed 
by noting that the last piece comprising the products of two /2’s and including an inhnite 
summation are nonsingular in e and can be neglected away from the crossover region. We 
then hnd at zero external momenta and special external quasi-momenta i the following result 


:xxx: 


(iv= + 6iV + 20), 
= 3 <t- - - 




+ 2eC<'>(!,r) 


(87) 


A similar observation permits our computation of the nontrivial two-loop diagram since 
the integrals I 4 do not contribute to the singular terms of this graph. At zero external 
momenta and nonvanishing external quasi-momenta set to the value i in all external legs, it 
turns out to be given by: 



,J5Af + 22) 

= OCT--- U 

27 



^ + 2eC*'’('i,f))}. 


( 88 ) 


In normalization conditions, we shall need the derivative of the two-point vertex part 
diagrams computed at zero external momenta in this massive framework. The “sunsetting” 
two-loop diagram involves the integral ^ _ o,i,r) and we are interested in 

^ d 

the object together with the other tilded integrals . The three-loop 

graph has a similar structure in terms of the integrals The 
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utilization of the combinations previously prescribed as well as the result of the integrals in 
the Appendix A imply that those diagrams can be cast in the form 


d 

d 


rO-. 


fc2=o 246 

+ 2){N+ 8) 




(89a) 


i-| + yir<'>(";f)]. (89b) 


fc2=o 162e2 

Some useful dehnitions similar to those occurring in PBC and ABC boundary conditions 
are: 




S, 


d 

j=-oo 


+ kj + i;r) 
[(g)2 + f2^-2 + l]/3 ' 




9 fS h I 


Fa,l{k,j,i;r) = ^f J (f 

.-r)- 

'’a KhJ^')— Qj^2 


,rf_i Fi^\q +k,r,r) 


[(g)2 + r2j2 + l]/3’ 


k=0' 


(90a) 

(90b) 

(90c) 

(90d) 


When expressed in terms of these dehnitions, the above hnite size correction to the higher 
loop two-point graphs represented by the amount can be written as 


= ---2 


1 rl 



dxdyil — y)ln 


0 Jo 


y{l-y)i +{l-y)r + 


yr 


-2 


1 rl 


+2 

+Qf 



0 Jo 


1 rl 



dxdyil - y)fi ( iy, ^2/(1 - y)F + {1 - y)r 2 + 

1 


yr 


-2 


0 Jo 
1 rl 


dxdyil-y)[r i {1 - y) + y 


a:(l — x) 


a;(l — x) 

-11+1' 


a;(l — x) 

-4F;,(v,f) 


+T / / dxdyil — y)i Cry + y 


0 JO 


x(l — x) 


1+1 


-12[7'f'(i,0if) + T7'^>(0,2if) 


(91) 


Finally let us present the solution of the vertex graphs at a suitable symmetry point 
(zero external momenta). These diagrams can be written in terms of those integrals of the 
four-point function. For instance, the one-loop diagram, according to the rule outlined in 
the diagrammatic Eq. (l38ll can be written as 

>50 


= T 


SP 


+ 


3p-^r (A^ + 2) 
18e 


t t t 


1 - - + 


• (92) 
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We are also interested in the solution for the two-loop diagrams. Using the same rule, we 
hrst obtain for the trivial two-loop diagram 


s: 


= T 


SP 


gw>y;=0 (^ j=2i 


(A^ + 2)2 

2 108e2 


[l-e + 2eC("nhr)].(93) 


The nontrivial two-loop, on the other hand is given by: 


(, 

A 


X 

, + / 

X 


V 


^ J SP c 


c 


f-. Z ( Z€ Z 


A glance in the results obtained from the graphs analyzed so far demonstrates the non¬ 
trivial character of the hnite-size correction terms for DBC and NBC in comparison to 
the BBC and ABC counterparts. The nondiagonal terms produce singular contributions 
at two-loop order in all primitively divergent vertex functions (and three-loop order in the 
two-point function) which modify the correction terms from the periodic and antiperiodic 
cases. 

With these data we can proceed to compute the critical exponents at least up to two-loop 
level. We carry out this task in the following section. 


IV. DBC AND NBC CRITICAL EXPONENTS IN THE MASSIVE FINITE-SIZE 
REGIME 


In order to compute the hxed point and hnally the critical exponents in a parallel plate 
geometries with one direction of hnite extent in the apropriate regime, we write the dimen¬ 
sionless bare coupling constant in terms of the renormalized one as = u(l+a!{^ u^). 

Although we are going to show that the hxed point is boundary condition dependent, we shall 
omit henceforth the explicit dependence on r in both bare and renormalized dimensionless 
coupling constants. The renormalization functions are written as = 1 -|- 
and = 1 + c^i'^u + 

The bare vertex parts which can be renormalized multiplicatively at zero external mo¬ 
menta and nonvanishing external quasi-momenta i can be written in the form 


af(2) 


dk"^ 

p(4) 


= 1 


M .,2 


b'PK + b^ 

\C) 


p(2,l) 


fc2=0 

= + (^2 
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SP 


(1 U _|_ ^(2 T )^„.2 
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OJ) 


_— —— C*i wq -l- {C 2 ’ C^ ' )u\ 


r^)ui 


(95a) 

(95b) 

(95c) 
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In the previous equations, the following identifications can be made: 


h{r) ^ 

2 6 

5(r) ^ ^ A 
^ 4 dk^ 


.-e. 

^ 2diA^J\i 


iN + 2) 


fc 2 =o 144e 

{N + 2){N + 8 ) 


1 - ^ + eW^^\i,r) 


648e2 

1 - I + 


ir>^C;xx>^ 


k ^=0 

{N + 8 ) 

6 e 

,• (iV2 + 6A^ + 20) r 1 


1-j + 


i 


(2 r)^ 


4(t 

cr 


36 


(l-e + 2eC^^Hh^))|> 


_ (5A^ + 22 ) / 

“ 18e2 I 


1 - I + 2eC^^Hb^' 


(7r) = 3p-(^ 


SP 




2 e 


SP 


SP 


3r (A^ + 2) 

1 6 e 
_ 3r (A^ + 2)2 
” ^ 36e2 

3t {N + 2) 
“ 1 12 e 2 


1 - - + 


[1 - e + 2 eC^'^^(f; r)] 


1 - - + 2eC^^^(i;h) 


(96a) 

(96b) 

(96c) 

(96d) 

(96e) 

(96f) 

(96g) 

(96h) 


Next, we apply the normalization conditions Eqs. fITOll . Using the dehnition of the renor¬ 
malized vertex parts discussed before, we can determine the normalization functions by 
imposing hniteness of the renormalized quantities. We start by replacing the dimensionless 
bare coupling constant in the diagrammatic expansion of the two- and four-point function. 
Use the dehnition from into Eq. flTObp in order to determine Replace the dimen¬ 
sionless bare coupling constant by the renormalized one at 0 {uq) since the correction will 
be O('U^) and can be neglected. We then End 



(AT+ 2 ) 
144e 


- + eW^^HiA) 
4 ' ' / 


(97) 


Using the above expression in the dehnition of the renormalized four-point function is suffi¬ 
cient to obtain 


(r) _ {N + 8 ) 


a = 


(t") 

02 = 


6 e 


’(Ar + 8 )' 

2 

6 e 



1 - - + eC^^)(f;r) 


1 -I- 2eQ^'^\i] r) 


{2N^ + 41Ar + 170) 


72e 


fr) 

Using these data, we can compute 63 ^ which can be shown to be given by 


6 ^ = 


(A^ + 2)(A^ + 8 ) 
1296e2 


1 - ^ + 3eC^^\i;k) 


(98a) 

(98b) 


(99) 
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The renormalization function which dehnes the renormalized composite held can be deter¬ 
mined in a very similar manner using the above expressions for and the diagram¬ 

matic expansion for that vertex part. We then hnd 


(X) _ (iV + 2) 


c ^ 


6e 


1 - o + 


(100a) 


(iV + 2)(iV + 5) (iV + 2)(2iV + 13) {N^ + 7N + 10) 

2 “ » 72i + 18i ^ 


The Wilson functions dehned in Eqs. fl75D can be written in terms of these coefficients as 


= —eu[l — (101a) 

= -eu[2b^2^ + (36^^^ - (101b) 

7^2^) = eu[c^i^ + ( 2 c^ 2 ^ - (4^V - (101c) 


In order to evaluate the critical exponents, we need the hxed point of the dimensionless 
coupling constant, which is determined by the condition /9*'^^(ttoo) = 0. It is a simple task 
to demonstrate that it is given by the following expression: 


^oo 


■ 6e ■ 

-ivTs- 


1 + 


[(9A^ + 42) 
- (A^ + 8)2 



( 102 ) 


The exponent rj is identihed with 'ylp\uoo) which implies the following three-loop result 



iN + 2) 6(3iV + 14) 

2(A^ + 8)2 I ^ (iV + 8)2 



(103) 


The other Wilson function at the hxed point is given by 7 ^ 2 ^ (woo) = [nXs) 
exponent u is related to the exponent rj and the last expression through = 2—("Uoo)—7, 


(A^+2) 


(6Ar+18) 

(W+8)2 ^ 

-, 7 ) I 


The 


which yields 


_ _ 1 , (iV + 2) ^ ^ (A^ + 2)(A^2^23A^ + 60)^2 
^ “ 2 ^ 4(A^ + 8)^ ^ 8(A^ + 8)3 ^ 

These universal results are independent of the boundary conditions and reproduce the 
bulk critical exponents. Despite the hxed point and the Wilson functions carry residual 
nonuniversal information from the normalization conditions employed encoded in the non¬ 
trivial correction function the cancellations take place exactly to eliminate all the 

hnite-size information resulting in the bulk critical exponents. A similar feature was shown 
to be valid for PBC and ABC. In the present case, even though the boundary conditions 
have not the same smooth character from PBC and ABC, we have just proven that the 
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previous conjecture made by Nemirovsky and Freed that the critical exponents are the same 
from the bulk theory also for DBC and NBC. 

In order to complete the analysis of the critical behavior from hnite-size systems conhned 
between parallel plates, let us discuss the situation from the viewpoint of masless helds, 
generalizing the picture just developed involving only massive fields. We shall use other 
set of normalization conditions and will show that the critical exponents agree with those 
obtained in the present section. 


V. MASSLESS FIELDS FORMULATION: DBC AND NBC IN THE EXPONEN¬ 
TIAL REPRESENTATION 


We have determined in previous sections the exact form of the several diagrams in per¬ 
turbation theory. Now we perform the substitution of the massive integrals by their zero 
mass analogues. Then, it is possible to define a simple set of normalization conditions which 
allow us to consider the same type of diagrams already considered in the massive framework. 

Let us get started by describing the normalization conditions. The restriction to the 
minimal number of diagrams is analogous to what was discussed in the massive theory. We 
dehne the quantity /i = 0,d) = /i = 0,d) — This means that the 

term does not need to be renormalized by the normalization functions, just like we 
did for the massive helds. (We careless kept the original vertex = 0,/i = 0 ,(t) in 

our treatment for ABC and BBC such that the similar term cr^r^ appeared therein. It 
has no effect for PBC, but produces a small deviation in the normalization as discussed 


in Ref. 


22l | for ABC. The correct form to get rid of this inconvenience is to dehne the 


object = 0,/i = 0,cr) = = 0,/i = 0 ,(t) — for ABC as we did above 

for DBC and NBC. This maneuver do not alter the normalization constants obtained in 
Ref. 22|, since the argument there implicitly took into account this feature). We choose 


a symmetry point with the following properties: i) since the theory now possesses infrared 
divergences, we have to renormalize the primitively divergent vertices at nonzero external 
momenta; ii) we also choose nonvanishing external quasi-momenta owing to the boundary 
conditions. If ki are the external momenta associated to the inhnite {d — l)-dimensional 
subspace and j is the moding attached to the external quasi-momentum characterizing the 
distance between the plates, the theory is renormalized at fixed external momentum scale 
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K and arbitrary nonvanishing moding i. The symmetric point is defined on the infinite 
subspace with the condition ki.kj = — 1), whereas all the external quasi-momenta 

mode of any primitively divergent one-particle irreducible (1-P/) vertex part is chosen in the 
value i (except for see below). 

The multiplicative renormalization can be successfully implemented through the following 
normalization conditions on the primitively divergent vertex functions: 


f^^\k = 0,i,9,0) = 0, 

(105a) 

dt^^\k = 

dk'^ k^=K‘^ ’ 

(105b) 

^^R{khii = i,g,^) = 3d^, 

SP 

(105c) 

= i,k 2 ,i 2 = i,Q,j,g,0) _= 3^. 

SP 2 

(105d) 


Note that the symmetry point is such that the insertion momentum is related to the other 
momenta in last equation through = {ki -|- k 2 )‘^- Second, recall from our discussion from 
the massive fields that the same rule applies for the vertex T^^’^\ki,ii,k 2 ,i 2 -,Q,j,g-, 0 )\^'- 
the tree-level vertex diagram corresponds to the combination 

3 3 3 

^ ~ ^ I(i=±(ii+i2) +^ l(i=o)) (106) 

where the signal of j is fixed. We attach to the external legs quasi-momentum i, such that 
the insertion quasi-momentum has two contributions: either j = 2i or j = 0. 

In order to get rid of the mass insertions, we can follow two different trends. The first one 
resembles the massive theory and can be formulated as follows. Start with a tree-level bare 
mass jif). Impose the condition T^‘^\k = 0,jUo,j,\) = Mo- Invert this equation to obtain 
ho(ho) just as we done in the massive case. This go could be interpreted as the shift in the 
bulk critical temperature. Replace this back into the diagrammatic expression of the vertex 
part = 0,/io,j, A) = This eliminates all tadpole diagrams. Next, express go{Mo) 

up to 0(Ao) and Taylor expand the remaining vertex part around k'^ = k^. The net effect is 
to eliminate all mass insertions in all primitively divergent vertex parts. Finally set go = 0. 

Although this argument makes perfect sense from the point of view of statistical me¬ 
chanics where the mass is identified with the reduced temperature (distance to the critical 
temperature), the reader working on quantum field theory might feel uneasy by starting 
with a nonvanishing bare mass, defining a three-loop bare mass and set it to zero after¬ 
wards, since it resembles a dynamical mass generation induced perturbatively. It is then 
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worthwhile to develop the most traditional argument in which the masss is zero in all orders 
in perturbation theory as we will describe now. 

Let us start directly with po = 0, which implies = 0,/io = 0,j, Aq) = 0. Conse¬ 

quently, one hnds the diagrammatic expression up to three-loop order: 




(107) 


In the right hand side are all tadpole diagrams up to this order. Since the tadpoles do not 
depend on the external momenta, we can replace this expressios directly in the diagrammatic 
expression in order to get rid of them. The right hand side diagrams computed at zero 
external momenta survive. The mass insertions cancel most easily: at one-loop order the 
condition = 0,/io = 0, Aq) = 0 = ^ implies that the contributions with 

tadpole insertions appearing in all primitively divergent vertex parts are zero identically up 
to the perturbative order considered. 

Plugging together all these facts into the bare vertex part T^‘^\k,j, po = 0, a) (for further 
details, the reader is advised to consult Ref. Q]), we end up with 


f(2)(A:,po = 0,Ao) = A:' 


6 




.(108) 


The situation here is completely similar to the massive case, since the number of diagrams 
are the same, only the renormalization symmetry point is chosen differently due to the 
infrared divergences occurring in the massless theory. We will keep the same normalization 
point for the external quasi-momentum modes i and a nonvanishing external momenta. 

An arbitrary vertex part including composed operators which is multiplicatively renor- 
malizable in a massless theory is dehned by the following expression 


T^i’'^\pn,tn,Qn',Jn', 9 , 0 ) = (Zj, ^, Aq, 0, A) , (109) 

where the quantities are the renormalization functions which make the above 

dehned renormalized vertex automatically hnite (even when the cutoff A goes to inhnity). 
As discussed previously our perturbative analysis will utilize only dimensional regularization 
and forget about the cuttof from now on, since formally the divergences of the integrals 
manifest themselves as dimensional poles represented by negative powers of the parameter 
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e = A — d. In order to prevent confnsion, in the above expresssion the identihcation = 


shonld be taken into acconnt. 


The external momentnm scale where the theory is renormalized induces a flux in space 
parameter such that the renormalized vertex parts satisfy the renormalization group equation 

~ ^ = 0 , ( 110 ) 


where /^^(n) = l'i\u) = and = /?(") • 


The 


dlnZ^A^ 


combinations zj^^ and 7^2^ = — g;;^) will be important in what follows. 

The renormalized theory is dehned at hxed L. The reader is advised to consult Refs. QQ- 


We turn our attention to the dimensional crossover regime in the massless theory. We 
will analyze the one-loop four point contribution and investigate its properties in the limits 
L —)■ CX3 and L —?• 0. The latter should indicate how small L can be without affecting the 
validity of the e-expansion. 

The graph required involves the combination of two integrals which follow directly from 
our discussion in the massive case, namely 


n 

l2{k,i,a) = J 


l = — CO ' 


h{k,i,j,a) = a 


[q‘^ + d^P] [{q + ky + d^{l + iY] 

dd-iq 


(111a) 


(111b) 


-|- a‘^P][{q + ky + ct'^P] 

We begin with the calculation of the the integral l 2 {k,i,d). Utilize a Feynman parameter 
X. Next integrate over the momentum using Eq. flA2j) from Appendix A. We are left with 

n 

a;(l — x){k‘^ + d'^P) 


ss,_ir(7i)r(2 - fcii) ^ 

Hk,i,a) = -^^— Y 

l = — OC> ' 


dx 


cr^(Z -I- ixY 


(d-l) 


-2 


( 112 ) 


Next perform the summation using the generalized thermal function as practiced before. 
Transforming Sd-i into Sd as previously prescribed, dividing the result by Sd and performing 
the e-expansion of the Gamma functions, we hnd 


a 


l2{ky,a) = -[l + -)[a^P + kY ^-h/ dxfi:e[ix,\x{l-x) 




,(113) 


In the massless theory, we factor out the external momentum scale. We dehne Fo,{k, i] a) = 

, factorize the external momenta by dehning the 


a 


-2a 


Jg dxfl_^_^ lix, W x(l — x) 


*2 + 11 
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quantity r = ^ and identify /2(«^,q5‘) = l 2 {K,i,f) (with similar identification for Fa) after 
this, which is the external momenta scale we need at the symmetry point. By expanding in 
e gives the following result: 





|/n(l + + ^Fo{K,r,f) 


(114) 


Following a similar trend it is a simple task to compute the nondiagonal integral whose 
solution can be written as 

rl 




rn 


dx[x{l — x) + F{Fx + j^{l — x))] 2. 


(115) 


'0 


In order to get the diagram computed at the symmetry point, we choose the external mo¬ 
mentum and quasi-momentum at hxed nonvanishing values. We have already discuss the 
combinations of last integrals entering in the process of computing the required diagram. 
After using standard manipulations (see for instance Ref. {mD, implies the outcome 


:o; 




1 + f + eC^^\r,f) 


= -Un{l + -F iFo(K,2z,f) -F iFo(«,0,f) 


-1-^(1 -I- 2r) arcsin 


'KVT 

3 


Vl+4r 


(116a) 


(116b) 


The hnite-size contribution function above has the same general structure as that from the 
massive theory, as it is going to be shown in a moment. Recall that the above coefficient 
a appears in the dehnition of the four-point renormalized vertex function. Just like in the 
massive approach, this factor will multiply the diagrams of this vertex function in arbitrary 
loop order and can be factored out from our asymptotic analysis. 

Before focusing on the asymptotical values for the correction function, it is important 
to mention that the would be scaling variable in the massless theory j goes to zero for 
hnite values of L as (^ = oo in that case. So the cases that could be assessed here are given 
respectively by the regions j < 1. The limit L —)■ cxd corresponds to j ~ 1, whereas small 
values for the distance between the boundary plates are characterized by the limit ^ 0. 

In the latter, the question again is how small L can be such that the e-expansion is still 
valid. 

We commence with the large L limit, which is the same as f —)■ 0. Since the linear terms 
in f goes to zero trivially, we restrict ourselves to the evaluation of the terms Fq{k, 2 i, r) and 
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Fo(/^)0,f). Using the definition along with the Boschi-Farina representation in terms of a 
summation of modified Bessel functions and taking that limit we find 

CXD /»! 


2i,f) = 4/ (ia;cos(47r?7-ix)A'o(27rn^ x{l — x)[r ^ + 4z^]) —»■ 4/ dxcos{4:TTnix) 

n=l n=l -^0 

pi 

X Kf){2nn^/x{l — x)r~^) < 4 / da;iiro(27rn\/x(l — x)r~^), (117) 


JO 


n=l 


where the latter is precisely Fo(K, 0 ,f), and we will show henceforth its convergence in 
that limit. Using the asymptotic values of the functions for large values of x, namely, 
^oiz) = + 0 ( 1 )), we have 


Fo(fi:,0,r) = 


1 oo ^ /-I 

2 f 2 y / 

^ Jo 


dx[x{l — x)] ^exp{—nB), 


(118) 


where B = 2'k^Jx{ 1 — x)f If we interchange the summation and the integral, we have to 
compute the summation J2n=i n-hxp{-nB) < exp{-nB) = ^ ^ < 

B~^. Therefore, 


r( 


1\2 


Fo(k, 0, f) < f2 —t 0. 


r(|) 


(119) 


Last equation is the rigorous proof that the correction goes to zero and one recovers the 
bulk result. 

Let us take the opposite limit f —)■ oo. We split the several pieces from which 
is made of and take this limit on each of them. We can easily obtain in that limit the 


asymptotic behaviors: —l/n(l + 4f^z^) —Mn{ri); ^(1 + 2 r) arcsin 
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\/l+4f2i2 




(r+2) 
3i ■ 


We are left with the nontrivial task of determining the terms proportional to Fo{k, 2i, r) and 
Fo{k, 0, f). Using Eq. flMD . it is actually simple to prove that \Fo{k, 0, f) —Mnr + ^f. 


Performing the identifications x = 2n and t = 


2ix 


Eq. flHTl) leads to the following expression 
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In this expression we can set r~^ —)■ 0. Although the dimensional crossover regime has to do 
with the behavior of large f (small values of L) we choose to keep the explicit dependence 
on i for the asymptotic behavior and neglect the other constants. The second and third 
terms are trivial in this limit, and yield the contributions 2ln2i and j, respectively. The 
two last terms involving the inhnite summation and the parametric integral are identical 
to those already worked out in the massive theory. In this limit they produce the quantity 

OD 

2 ^ i^"(C(2n+1). Dehning again the regularized correction for the massless theory : ([i] f) : = 

n=l 

^ OD 

C(i; h) — I X] 1) which is identical to that for small values of L in the massive theory. 

n=l 

Altogether, the regularized finite-size correction for small values of L (r — )■ cxd) in the massless 
case reads 


Afr)/- , '^r{l + T) (2r-h5) 

: Q ’{i]r) —Inr -\ -- - + 


( 121 ) 


3 6 i 

It is simple to see that the hniteness correction is compatible in massive and massless cases. 
Notice the absence of terms proportional to Ini which got cancelled along the way in both 
massless and massive theories. Since the presence of this sort of term in conventional field 
theory defined for inhnite systems implies that all lenght scales are coupled, the absence of 
them in the mode of external quasi-momentum i shows that there is nothing fundamental 
on the sector of the hnite-size correction which depends on i. In other words the layered 
geometry do not themselves interact strongly up to all length scales, from one plate to the 
total number of them, since this information is encoded in the logarithm which is absent in 
the correction. The apperance of the relevant lenght scale L marks a profound difference with 
the bulk critical behavior even though the expoents are the same. It is a good explanation 
why surface effects are not important in this regime. 

If we choose the Neumann boundary condition r = 1, we hnd the asymptotic behavior 

27rf 


: r) —Inr -|- 


3 ’ 


( 122 ) 


and the leading divergence on r is linear as happens with PBC. This is in agreement with 
the massive case. 

On the other hand, if we take Dirichlet boundary condition we have instead 


: Q'^\i', r) —Inr, 


(123) 


and the divergence is purely logarithmic as that for ABC, which is just compatible with 
the massive case as well. If L is too small, those singularities take over the poles in e of 
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the one-loop four-point function therefore invalidating the e-expansion results. This is the 
meaning of the dimensional crossover regime, where hnite-size correction effects get bigger 
than the leading singularities in e as already discussed for PBC and ABC in Ref. 22[ |. 


The reader might have noticed that, except for the ^-dependent terms appearing in the 
hnite-size correction function in the L —?■ 0 limit, the rigorous arguments used to determine 
the correction function for PBC and ABC carries over for DBC and NBC. Since we are 
neglecting the dimensional crossover regime, that is the main reason we refrain from writing 
another appendix with the derivation of the results for the massless integrals. 

With the resources furnished in the present paper in conjumination with those details 
worked out in Ref. 2^ for ABC and PBC, the other two-loop diagrams of the four-point 


function at the symmetry point can be shown to be given by 


:xxx: 


= 3a 


{N^ + 6N + 20) _2j 1 
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Similarly, the derivative of the two- and three-loop diagrams (with respect to at the 
symmetry point) of the two-point vertex part can be computed similarly. From our discussion 
so far, it is not difficult to hnd out the following (divergent) expressions 


d 

d 




k-^^{N + 2 ) 


fc2=0 


24e 
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l + — - 2 €W^^\i,r) 
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(125a) 

(125b) 


’fc2=o 162e2 

Here the hnite-size correction for the two-point function as well as the quantities which 
dehne it can be written in the form 


H/T)(^; f) = -f 2F'o{k, i] r) — Fo(«^, h r) — 3f[Ho{K, i, 0; f) 
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.The other amounts appearing in the equation for the 


where 

correction are dehned by 


fc2=K2 








[g^ + a'^FY 


'^a 1^5 ’) — Qj^2 


fc2=^2 * 


(127a) 

(127b) 

(127c) 


We have all elements to compute the critical exponents using the massless helds frame¬ 
work. As we are going to show in the next section, all dependence on the hnite-size cor¬ 
rections for the several vertex parts disappears in the Wilson functions at the nontrivial 
attractive infrared hxed point. The main consequence of this is that in the hnite-size regime 
the critical exponents are boundary condition independent. 


VI. CRITICAL EXPONENTS IN THE MASSLESS SECTOR 

We are going to be rather brief in the present section. We shall focus solely on the 
computation of critical exponents using normalization conditions, albeit the results in the 
present work should be worthwhile in the evaluation of critical exponents using minimal 
subtraction. The algorithm was already set in the massive case. All we have to do is to 
replace the massive diagrams by its massless counterpart. It is important to realize that in 
the perturbative expansion of the primitively bare vertex parts, the bare coupling constant 
is dehned in terms of the bare dimensionless coupling constant through Aq = This 

implies that for arbitrary hxed n, the perturbative expansion consisting of the multiplication 
of suitable powers of Aq by the diagram eliminate all the k prefactors appearing in each 
diagram as just discussed in the previous section. 

Firstly, use the expansion of uIY\ and in powers of u in exactly the same form 
as we did before in the massive case. Secondly, employ the analogue of Eq. fl95ll recalling 
that now the symmetry point occurs for nonzero momentum and the diagrams correspond 
to massless helds. We refrain from giving the massless version of Eq. ([96]), rather we 
implicitly plug them in the dehnition of the primitively divergent massless vertex parts. We 
then obtain the following solution for the renormalization functions (i.e., and Z^2) 
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coefficients 
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^ + 2)(^ + 5) 2N^- + lliV + 14 + 7iV + 10 

2 36e2 72e 18e ^ v ; v ; 


We are going to use this information in the calculation of the Wilson functions. We just 
have to plug these coefficients into Eqs. fllOip . The nontrivial infrared fixed point value of 
the coupling constant u* is obtained from the condition j3^'^\u*) = 0. It is given by 


u = 


6e 


W + 8, 


1 + 


(9W + 42) 


o -C 


W| 


i\r] 


(129) 


(N + 8)2 2 

Upon substitution of the hxed point into 7 ,j)(m) one obtains the exponent rj (the anomalous 
dimension of the field) up to three-loop order, namely 


_ f *\ ( 2 ^ + 2) 2 

T] = 7^(m ) = :rTT^-rr^e 


1 + e 


6(3W + 14) 1 

(W + 8)2 “ 4 


2(W + 8)2 

Finally, using the relation = 2 — 7 ^^(w*) — t], leads us to the result 


(130) 


u = 


{N + 2) (A^ + 2)(W + 23W + 60)^2 


(131) 


2 ' 4(A^ + 8 ) 8 (W + 8)3 

In spite of being quite different as far as the e-expansion of the diagrams are concerned, 
the massless and massive settings yield the same critical exponents as expected from the 
universality hypothesis. Ultraviolet and infrared analysis are completely equivalent in the 
evaluation of critical exponents, which reproduce those from the bulk as in PBC and ABC 
conditions. 

What is really appealing in this new formulation inspired in the Nemirovsky and Freed 
construction for massive fields is that no surface fields are required to implement successfully 
the hnite-size renormalization program for DBC and PBC. 
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VII. DISCUSSION 


Formerly, there was a certain misunderstanding by some authors concerning the phe¬ 
nomenological hnite-size scaling hypothesis. It was thought that the variable j governed 
the approach to bulk criticality: whenever j —)■ cxd the e-expansion results shown for the 
above diagram could be trusted, whereas j < 1 would be the regime in which perturbation 
theory could no longer be reliable. In the modern approach to hnite-size criticality proposed 
in the present work, this scaling variable looses its meaning as discussed in the ABC and 
PBC cases. In fact, we are going to show next that even in the massive case for DBC 
and NBC, the regions j ~ 1 and j < 1 are also available to scrutiny within the present 
technique. 

How small j (or how large r) can be without invalidating the e-expansion results? Setting 
yU = the variable can be rewritten as r = so that the limit L —)■ 0 corresponds to 
r —)■ cxD. Before we start our discussion, identify the dimensional single pole (oc e“^) with 
the = ln{A^) in a theory regularized with a cutoff, which will be helpful to us in what 

follows. 

The dimensional crossover condition for NBC implies that the linear term in r is greater 
or equal to that coming from the pole in e. (We can extract further information by identifying 
the ultraviolet cutoff with the inverse of the lattice constant a as A ~ -.) In other words. 


L < . 

3ln{^) 

Away from the critical temperature ^ is hnite and typically a multiple of the lattice 
parameter. Take ^ = 10a. In that case, L < 9.1a. There is a narrow region 9.1a < L < 10a 
where j < 1 and the perturbative expansion is still valid. The number of parallel plates 
is n = ^ -|- 1, where n is an integer. Then, for n < 10.1 the e-expansion results break 
down. Since n has to be integer, this condition implies that for n < 11 the perturbative 
expansion is invalid. However, it is valid (at least formally) for n > 12. Hereafter, anytime 
we mention the region of validity of the e-expansion for L in a certain range, it is implicit 
that this variation interval on L obeys the criterion j < 1- Now take ^ = 100a, which implies 
L < 45.5a, n < 47. There is a wider region 45.5a < L < 100a which does not invalidate 
perturbation theory. For ^ = lO^a, the breaking conditions are L < 303.2a (or n < 305) and 
the window of validity of perturbation theory is contained in the range 303.2a < L < lO^a. 
By allowing ^ = lO^a leads to the condition(s) L < 2274a(, n < 2275) and a wider range 
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(2274a < L < 10000a) exists which does not invalidate perturbation theory, and so on. 


In PBC from Ref. 


22|. the linear term was equal to where r = Using exactly the 


same argument as above in NBC, the dimensional crossover condition reads now L < , . 

- MU 

When ^ = 10a implies L < 13.64a and n < 15. For that value of the correlation length, 
perturbation theory is only valid for y > 1- The choice ^ = 100a turns out to produce the 
breaking conditions L < 68.21a and n < 70, with a region 68.21a < L < 100a where j < 1 
does not invalidate e-expansion results. Taking ^ = lO^a, perturbation theory is invalid for 
L < 454.8a or n < 456. The window of validity is 454.8a < L < 1000a. Finally, the choice 
^ = 10‘^a leads to the dimensional crossover conditions L < 3410.9a and n < 3412, with 
a wider region of validity, namely, 3410.9a < L < 10000a. Comparing with NBC results, 
the bulk critical behavior requires a smaller number of plates for NBC than in PBC and 
a thinner him geometry. Since fewer plates are required in NBC to keep the bulk critical 
behavior, from the energetic viewpoint it is easier to provoke a bulk-surface transition in 
NBC than in PBC. Therefore, our approach explains why nature chooses NBC as the 
prototype for the bulk-surface transition, since this information is encoded in the hnite-size 
correction to the bulk behavior. 

We turn our attention to the dimensional crossover regime for DBC. The logarithmic 
divergence in r —>■ oo in DBC can be understood similarly: regularizing the theory with a 


cutoff A such that 


In 


The correction will become as big as the dimensional pole 
whenever L < vra ~ 3.14a, which is independent of the bulk correlation length U Using the 
same line of thought for ABC yields L < 2na ~ 6.28a for the collapse of the bulk critical 
behavior description in perturbation theory. Both results are independent of the (fixed, 
finite) bulk correlation length U For different fixed bulk correlation lengths one needs n < 5 
(n < 8 )for DBC {ABC) by demanding integer values for n. For n > 6 (n > 9) the 
bulk critical behavior dominates the finite-size corrections and there is a real experimental 
possibility of construction of nanodevices for DBC {ABC). 

For instance, by slicing a thick material presenting bulk critical behavior in thin films 
satisfying the above conditions, one could vary the temperature to different values from the 
bulk critical temperature (variation of the correlation length) and use neutron scattering 
experiments to see whether the bulk behavior varies for large enough values of n. If it does 
not, the boundary conditions on the limiting plates are either DBC or ABC and we have 
a practical way to determine the boundary conditions in a certain material. If it does, the 
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boundary condition should be either NBC or BBC and the number n should be much 
bigger than in the previous situation. 

The description of the dimensional crossover regime in the massless case is also worth¬ 
while. It is important to mention that hnite values of L already implies 0 since .^ = cxd in 
the massless (critical temperature) formulation. The important quantity here is the external 
momenta scale which generates the renormalization flow, although it has no obvious physical 
interpretation like the connection of the mass with the correlation length. Nevertheless, we 
can vary it and see what happens with the correction. 


For NBC, the linear term in the massless case is where f = - (a = f) as can be 
verihed from Eq. (11211) . What is left here is to perform the variation of k in units of the 
inverse multiple of the lattice parameter. Following the same trend we hnd the critical value 
For instance, we can set = 10a, = 100a, = lO^a, and = lO'^a. 




L < - 

3«:/n ( — ) 

Obviously, the picture does not change. All the critical values of L below which the system 
undergoes the dimensional crossover gets a factor of tt with respect to those from the massive 
theory. For instance, kT^ = 10a implies that L < 28.6 (n < 30); for = 100a, one has 
L < 142.9 (n < 144); when = lO^a, we obtain L < 952.5 (n < 954); etc. 

For BBC, a comparison involves the correction to the four-point function which is given 
by the second, third, fourth and hfth terms from Eq. (61) from Ref. 22|] by taking r = 0. 
Of course, we have to reconstruct that expression (taken at k = 1) by multiplying all terms 
proportional to L by k. Then, the second term gives the logarithmic divergence on r = ^ 
(a = ^) while the third gives the linear divergence on r and the other terms are not 
important in what follows. Altogether, the aforementioned correction for small values of L 
has the linear divergence as the dominant term which invalidates perturbation expansion 


whenever L < 


It is not necessary to go into the details, to see that several values 


Kln( — ) 

of K in the same range as before produce critical values of L which are bigger by a factor of 
TT in comparison with the analogous situation in the massive theory. The reader can check 
that the same feature takes place to NBC when compared to the massless case. This should 
be not surprising since the fluctuations are enhanced in the massless regime which requires 
thicker slabs to guarantee the bulk criticality. 

Perhaps the most interesting result regards DBC and ABC which have the same critical 
value for L in both massive and massless regimes. This is so because in the massless and 
massive regimes the logarithm term for the variables in both situations does not change its 
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coefficients. The full unveiling of this invariance remains to be investigated. 

In closing this discussion, if these inequalities are not satisfied, the new situation corre¬ 
spond to new critical behavior not belonging to the usual bulk universality class, since the 
smallness of this variable disturbs the system in an uncontrolable fashion driving it to a 
“dimensional crossover”. Note, however, that the regions ^ ~ 1 and y < 1 are away from 
the dimensional crossover regime since they are not too close to zero. They can be safely 
described within the present modern hnite-size scaling approach. 


VIII. CONCLUDING REMARKS 


In this work we have developed a momentum space method to calculate critical exponents 
for a critical system whose order parameter is dehned on a layered geometry with Dirichlet 
and Neumann boundary conditions imposed on the limiting (parallel) plates at z = 0, L 
with z characterizing the finite-size direction perpendicular to the {d — l)-dimensional (hy- 
per)planes of infinite extent along all its linear dimensions. The main features introduced 
are: i) the tensors needed to construct arbitrary loop diagrams for IP I vertex parts through 
the utilization of the exponential representation for the basis eigenfunctions with or without 
insertion of composite operators which can be renormalized multiplicatively; ii) the repre¬ 
sentation of the finite-size correction written in terms of a sum of Bessel functions instead 
of that involving integrals without elementary primitive; Hi) massless framework for the 
computation of critical indices. 


There are many similarities with the PBC and ABC treated recently 22|. For instance, 
the unification of the results for the boundary conditions, except that for DBC this only 
occurs as long as the external quasi-momentum are set to nonzero values. (For PBC and 
ABC this unihcation was obtained using vanishing external quasi-momentum.) This re¬ 
striction causes no loss of generality, since it is possible to show that NBC can also be 
formulated with zero external quasi-momenta and the results do not change. We leave this 
topic for a future work. The region of validity of the e-expansion is consistent either in the 
massless or in the massive formalism. In that region hnite-size scaling is valid, the dominant 
critical behavior is the one from the bulk with respect to the critical exponents even though 
the hnite-size corrections might appear in other universal quantities, like amplitude ratios 


32, l33j. From the dimensional crossover regime viewpoint where the e-expansion results are 
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no longer applicable when L gets smaller and smaller, DBC and ABC diverge logarithmi¬ 
cally with L whereas NBC and BBC diverge like L~^. This conhrms a previous claim by 
Nemirovsky and Freed in their Green’s function framework for the massive case. 

When the e-expansion ceases to be valid, those hnite-size corrections become bigger than 
the poles in e (representing logarithmic divergences of the bulk theory). In that case, the 
correction for NBC oc L~^ in the limit L —)■ 0 dominates over the term InL (which is one 
piece of the total contribution). In the present framework, the onset of the dimensional 
crossover is directly connected to surface effects but not in the way previously imagined, 
i.e., the passage from (A^, d) to the (A^, d — 1) universality class: the interaction between the 
plates provoke a more complicated effect transliterated in those L-dependent terms. In NBC 
the dimensional crossover regime can be identihed with the bulk-surface transition region 
although the identihcation is not complete since a proper description of this region should 
probably include external surface fields as well [l2, 34) 35|. In other words, the dimensional 


crossover regime marks the transmutation of hnite-size contributions into surface effects. 

From the technical point of view, there are other similarities among [DBC, NBC) and 
{ABC, BBC). Our utilization of the basis functions in terms of the exponentials for {DBC, 
NBC) makes it possible to use results from {ABC, BBC) analysis like performing summa¬ 
tions in the range (—cxo, cxo), since the functional form of these results are invariant in spite 
of minor modihcations in the hnite size parameter r (r,f, etc.). Besides, the terms which 
have been discarded along the way for being more regular than those kept in all stages of 
the process revealed themselves identical in their functional form for all of these boundary 
conditions, with differences obviously in the massless and massive approaches. 

But there are also totally different aspects when comparing {DBC, NBC) and {ABC, 
BBC). The latter are boundary conditions whose terms do not break translation invari¬ 
ance. The former actually do have translation invariance breaking terms. Previously, the 
common belief about translation invariance breaking was necessarily attached to surface 
contributions. These could stem from two origins, namely, if either surface helds i) are not 
allowed to begin with and these contributions are corrections to hnite-size effects or ii) are 

n n 

permitted and the subject goes beyond the hnite-size problem itself [12|, [1^. According 
to our study here the description of hnite-size systems subject to DBC or NBC can be 
understood entirely out of helds within the volume between the plates without necessity 
of referring to surface helds. Plus, the breaking of translation invariance in our treatment 
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has nothing to do with surface effects, provided the critical system is kept away from the 
dimensional crossover regime. This is highly desirable in comparison with the deconhned 
situation: the approach to bulk criticality in the region 0 < ^ < cxo clearly leads to a novel 
modern hnite-size scaling paradigm not arising from externally imposed surface helds. The 
structure unveiled in the present work demonstrates that translation invariance breaking 
comes from the “nondiagonal terms”. These terms actually make it difficult to get rid of 
mass insertions in order to have a smaller number of diagrams in the computation of critical 
exponents, for instance. Our prescription to eliminate the mass insertion graphs without 
invoking external surface fields or any sort of fields other than that representing the confined 
bulk system might be worthy to take the subject to another level of understanding. 

One interesting aspect is to pursue other renormalization schemes like minimal sub¬ 
traction for massless and massive fields for the boundary conditions just presented. The 
discussion for massless fields subject to ABC and PBC already appeared in Ref. |22|. A 


thorough discussion for all boundary conditions utilizing massless and massive fields is still 
lacking. It would be nice to see whether setting the external quasi-momenta at fixed values 
should be sufficient to renormalize the massive theory 251. 


It is tempting to employ the machinery 


m-fold Lifshitz type competing systems 


36 


ust developed to treat the case of anisotropic 


37l |. A simple follow-up idea is to take the 


finite direction perpendicular to the m(-dimensional subspace) competing axes and investi¬ 
gate whether the critical exponents are affected by the introduction of this new ingredient. 
Amplitude ratios such as that from the susceptibility (see Ref. |38| for the m = 1 case) and 
specific heat [39| could also be tackled within the technique introduced here in order to figure 
out whether the finiteness modify them in comparison with those from the bulk. To extend 
this topic to its full generality, we could use it in the investigation of finite-size corrections 
of arbitrary anisotropic competing systems 2, Al\ and take the finite-size direction perpen¬ 
dicular to all types of competing axes inherent to the problem. The treatment of amplitude 
ratios for generic competing systems, for instance, those discussed in Refs. 42|, l43| poses no 
obstacle in principle and can be investigated too. 

There are many applications in the critical phenomena context which could be unraveled 
utilizing the results contained in our present results. For instance, further understanding 
of this topic would lead to the fabrication of new devices involving thin films of materials 
displaying bulk critical behavior. The importance of our work in guiding this enterprise to 
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experimentalists is the indication of how small L can be without spoiling the bulk critical 
properties, which is certainly important in the nanotechnology scenario. 

Since we have shown that the critical exponents for a noncompeting system in a parallel 
plate geometry with one direction of hnite extent is identical of those from the inhnite system 
we could try to hnd an exact perturbative analytical solution to the uniaxial (m = 1) 
Lifshitz critical exponents. In fact, if we allow the hnite-size to be along the competing 
axis the problematic quartic integral would get transformed into a summation which is 
easier to perform at least in principle. Taking this achievement for simple enough boundary 
conditions like PBC would automatically yield the critical exponents of the bulk systems 
for the uniaxial case. Can this be possible? The results just obtained will help to make sure 
whether this idea is feasible in the near future. 
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Appendix A: Massive integrals in dimensional regularization 


In the text we have omitted some steps concerning the one-loop integral l 2 {k, i, r) in order 
to get Eq. (175]) due to its similarity with PBC and ABC integrals already worked out in 


Ref. 


22l | . We commence by deriving explicitly Eq. fl75]) along with the nondiagonal integral 


J 2 (A;, i, d,/i). We shall plug the renormalized mass p in all diagrams to be discussed for 
the reasons discussed in the main text. 

Let us get started with Eq. fl26ap . We factor out the renormalized mass p and utilize 
utilize a Feynman parameter to obtain 


OO p 

= r/i"" J dx J 


d^-^p 


j=-oo 


X 


p^ + 2xkp -|- xk‘^ -|- C[{j -|- ixy + x(l — x)P] -|- 1 


(Al) 
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Dimensional regularization is usually expressed in the Feynman’s integrals evaluations 
through the identity (see Ref.j^) 


d'^q ir(f)r(a- 




(A2) 


J + 2A:.g +/i^)" 2 r(a) 

where Sd is the area of the d-dimensional unit sphere. Replacing this into the previous 
expression leads us to 




••I OD 


X dx ^ [a:(l — a:)(/c^ + + r^(j + + 1] 2 
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j=-oo 


The summation can be performed using the generalized thermal function identity 

00 

Da{a,b)= [(u + a)2 + 


24| 


- ^ 
r{a) 


n =—00 
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+ faici, b) 


where 


CXJ 1 

fa{a,b) = 4 cos(27rma) ^ 


m=l 


(A4) 


(AS) 


and K^{x) is the modihed Bessel function of the second kind. By identifying a{x) = ix, 
b{x) = + PP) x(l — x) + 1 and performing the continuation e = 4 — d, we hnd 




X 


r(-)[x(l - x){k‘^ + zV) + 1 ] 5 + + . (a, b) 


(A6) 


Transforming the argument of the T-function with the recipe y^r(^^)S'd_i = r(|)S'c; along 
with the expansion in e = 4 — d, we get to 

I^\k,i]d, ix) = Sdfx~^ -^1 — -') X f dx[x(l — x)(fc^ + z^r^) + 1]“2 

Lev 2 / Jo 

-1 

+-r“T^2 —J dx/i^zx, f~^\/x{l — x){k‘^ + PP) + 1 j . (A7) 

After dividing by Sd this object is precisely l 2 {k,i',r) from Eq. fITS)) (see also Eq. fl76)) b 
Setting k = 0 and neglecting 0(e) contributions leads to the following equation 


^ 2 ^\k = 0,i]r) = /i l — S — 1 I dx/rz[x(l — x)z^r^ + 1] + hFo(0, z; r) . (A8) 


If, e e 


2 2 
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We also need the integral l 2 {k,i, j; a, fi). Factorizing the mass and using a Feynman 
parameter in this integral from Eq. fl26bp one hnds 


= rfi " ^ / dx I d'^ 


OO pi 


j=-oo 

1 


X 


+ 2xkq + xk? + r‘^[{Px + (1 — x)p] + 1 


2 • 


(A9) 


Integrating over q as before, using y'^F(^^)S'd-i = F(|)S'd and performing the continuation 
e = 4 — d, we can rewrite this object as 

I^''\k,i,j-,d,fi) = [ dx[x{l - x)k^ + [xi^ + {1 - x)f]r^) + 1]~^ , (AlO) 

^ Jo 

that results in Eq. fl77|l after dividing the last expression by Sd- In particular, the combi¬ 
nation which appears to produce the one-loop diagram with all external quasi-momenta set 
to i is given by 




.(A^ + 8) 


9 


+4:l2{k,i,i;f) 


hik, 2i; r) + 2 l 2 {k, 0; r) 2 Tl 2 {k, 0, 0; r) -h 2 Tl 2 {k, 2i, 0; r) 

(All) 


Utilizing the results above mentioned evaluated at A: = 0 yields Eqs. fl79D . This completes 
our Erst task. 

Let us consider two- and three-loop integrals. We preclude the dimensional crossover 
region for too small values of L which spoils the e-expansion analysis. 

The easiest two-loop contribution for the four-point vertex function is given by the dia¬ 
gram which consists of the following combination of integrals 


ll{k, 2i; a, p) -F 2/3 (A:, 0 ; a, /i) -h ^Thik, 2i; a, /i) 


~ + 6A^ + 20) 

.XJUxi” ^ 27 

xl2{k, 0, 2i; d, /i) Arhik, 0; d, /i)/2(^, 0, 0; d, /e) -h 4J2(fc, i, f; d, ^)[l2{k, 2i] d, /r) 

CXD 

-h/2(fc,0;d,/i)] + 2 ^ [il{k,j,j + 2i;d,fx) + ii{k, j, j; d, fx) + rhik.jj + 2z;d,/i) 


j=-oo 


'^h(k,j + s,i + 'kcr,ij,) + -i,j + 


(A12) 


computed at zero external momenta. Since the contributions inside the summation in the 
above expression is regular (0(e°)), we can neglect them. It is then easy to show that 

1 




27 


Le^ 


1 — e -|- 2e( 


d)i 


i:r 


(A13) 
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The nontrivial two-loop diagram contributing to the four-point vertex function with ar¬ 
bitrary external momenta was derived in the main text and corresponds to Eq. fl5^ . The 
integral k 3 ,i, j, k, a, n) (where P = ki + /C 2 ) is regular and does not contribute to the 
singular part of the diagram. At the symmetric point it reads 

(5A^ + 22 ) 




27 


-a 


h{0, 2i, i) + 2 / 4 ( 0 , 0 , i) + 8 / 4 ( 0 , i, z, 0 ) 


-l-2r/4(0, 0, 0, i) -|- r/4(0, 0, 2i, i) -|- r/4(0,2z, 0, i) -|- /4(0, z, z, 2z) -|- 0{ii) (A14) 

The diagonal terms are composed by the integral / 4 ( 0 ,z, j) which can be written as 

■d-i 


jl = -QO 

When we factorize p we hnd 


d 


[q^ cr2(ji - z)2 + a'^jl + ^2j 


(A15) 




-2e 




Jl=-oo 


X 




[q,2 _|_ _ ^)2 _|_ lj^q2 _|_ j.2j2 _|_ 

(A16) 


where 


F±,2{0,i,j)= \ 




[q,2 _|_ _ j^^2 _j_ ]^j|'g2 _|_ ^2j2 _|_ 


(A17) 


Jl=-oo 

The singular terms in come from the high-momentum region of the momentum integra¬ 
tions, since they correspond to ultraviolet divergences. Focusing our attention in the last 

term, it will contribute to the singular part if and only if the object F^\q,ji — j) is pro- 

2 

portional to q^,p > 0 in the limit q —)■ 00 . We can discard it in the computation of the 
singularities if we can prove that it is proportional to q^ for p < 0 in the ultraviolet region. 
Indeed we can neglect this term as follows. 

First, write it in the form 


F±\q,j',r)=4:r / dx'Y^cos[27rmjx)] 


m=l 


TTUl 


r~^^/x{l — x){q^ -I- r‘^p) -|- 1 


X ^x{l — x){q^ + ^2^2^ _|_ j _ 

It follows trivially that 


(A18) 


= 4r" 

2 


dx 


nm 


— X)(g2 + f2j2'j _j_ I 


X K^i^nmr ^ y/x(l — x) (g^ -|- r‘^p) -|- 1 j > F^\q,j]r). 


(A19) 
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We analyze the latter henceforth. Since the integrand is symmetric aronnd x 
write 


nm 


ff\q,3;r) = Sf W ^ - w ^ , -2 -2^ , 1 

^ Jo 0- ^^^x{l - x){q^ + r^j^) + 1 

X (^7Tmf~^\/x{l — x){q‘^ + f'^p) + 


we can 


(A20) 


When g —)■ 00, we can take a small real parameter A << 1 with the property \q^ —)■ 00. We 
split the integral into two pieces: in the hrst one we integrate in the interval ( 0 , A) neglecting 
the term in the integrand, whereas in the interval (A, |) we replace the Bessel fnnction 
by its asymtoptic form. We then write 

00 

q—>-oo 2 


dx 


Vr-i-^x(g2 + r2j2) + 1 


X K^(2Timr ^\/x{q^ + r'^p) + 1 j + 


dx 


Tim 


X 


TT 


dvrmr ^y^x(l — x){q^ + + 1 


exp 


A '^r — x){q‘^ + r'^p) + 1- 

—27rmr“^y/x(l — x)(g^ + f'^p) + I'j .(A21) 


The second term goes to zero exponentially and can be disregarded. After replacing the 
change of variables y = 1 + x{q^ + in the integral, we hnd 


limfy\q,j]r) 

g^cs:) 2 


We also know that 



8r 


” m=l 


-^(7rm)2 / dyy 4i^|(27rmr ^y^). 


(A22) 


dxx 2(x —l)'^ ^iA',^(a\/x) = r(/z)2'^a (A 23 ) 

Considering the region ontside the dimensional crossover region where r < 00 implies the 
following asymptotic expression expression 



limfy\q,j]r) 

q^oo 2 


8 r^ e 



00 

^(7rm)^“^A'l_i(27rmf“^). 

m=l 


(A24) 


Since fy\q,j]r)(> F^\q,j-,r) is regnlar in e, in the nltraviolet region F^\q,j]r) and the 
2 2 2 

aforementioned integral involving it are both regnlar and we do not have to worry abont its 
contribntion to the singnlar part of the integral 14(0, z,j). 
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We can write the remaining expression in the form 




/ I oo p 

dx[x{l-x)]^ j (f-^q 


Jl=-oo 


1 


X ■ 


[g 2 + Pjl + 1 ] [(g 2 + ^ ' • 

We introdnce the Feynman parameter z to rewrite last eqnation as 

,1 oo „i 


[g 2 + r 2 (ji - z )2 + 1 ] 

(A25) 


4^^(0,hi) = ^'-(l - I) / dx[x{l-x)]^^ 


Jl=-oo 


d‘^ 


J [g2 + r2j2 ^ ^ ^ _ 2j^)^]2[(g2 ^ ^ y(^]^ ’ 

followed by another y to melt the two denominators into a single one. We then obtain; 


(A26) 




- 2 .r (2 + 5) 


r(i)£ 




Jl=-oo 


d^ ^q 


dyy{l-yY^ ^ 

(A27) 


J [g2 + r‘^jl + r2[i(z - 2ji)yz + (1 - y)j{j - 2ji) + y + 

Employing Eq. flA2ll to perform the momentnm integration as well as onr recipe to transform 
the nnit sphere area S'^-i into Sd, one can show that 


/r(o,2.i)=t^hiz|£(i±i)(i_£ 


dx[x{l — x)] ^ / dyy{l — y )2 


^-1 


2v^ r(|)£ 

OO ^ 

Y + (1 - y)f - ‘^jiiwz + (1 - y)j)] + y+ " • (^28) 

— <-V-\ ' ' 


Jl=-oo 

By evalnating the snmmation using the generalized thermal function upon the identihcations 
a 

factor Sd, that integral reads 


>y CVcllUcLLlli^ Lliti bUliilliclLlUli Ublli^ LliC clllZitiU LliCiliicll iUliULlUli UpUli ILLti lUt^liLliiUcLLlUiib 

= -iyz - {1 - y)j^ b = y Pyz{l - yz) + fy{l - y) + f-‘^ (y + and absorbing the 


iY\^S.3) = 


2 ^ 


i_£Nr( 2 -|) /•' 


-1 


2 ) r(f) 


X 


T{e)[i yz{l-yz)+j y{l-y) + r [y + 


— I dx[x{l — x)] 2’' I dyy{l — y)^2 ^ 
1 - ?/ 


'0 


'0 


a;(l — x) 


+fUe - (1 - ?/)i> t/ i‘^yz{l - yz) + fy{l - y) + r ‘^(y + 


1 - 1 / 
x(l — x) 


The integral over y presents a pole in y = 1. We then evaluate it at y = 1 j^, which is not 
only a lot easier but maintains the essential pole contribution of interest. Expanding in e 


.(A29) 
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and getting rid of 0(e°) terms, we finally obtain 

- I - e f dzln[r^i‘^z{l-z) + l]+eFo{0,r,f)]. 


(A30) 


We still have to compnte the nondiagonal contribntion coming from 14(0, It can be 

expressed in terms of the one-loop integral l 2 {k,i,a, n) as 


= a d q - 

J [q^ + + fi^][q^ + aH^ + 12^] 

Extracting the mass from the integrand we get to 


^4(0,hi, 0 = ^/^ / d ^q 


(A31) 


(A32) 


[g2 _j_ j.2j2 _j_ ]_j|'g2 _|_ j.2j2 _|_ 

Note that the integrand in the last expression is symmetric by the change i ^ j. This 
integral differs from 14(0, z,j) becanse i) the snmmation over ji is missing, ii) in the one- 
loop snbdiagram there appears only a external qnasi-momentnm I, Hi) each propagator 
in the “external bnbble” are attached to external qnasi-momentnm. The compntation is 
entirely analogous with the above computation for 14 ( 0 , i,j). Picking out only the singular 
term and putting aside the regular terms we are led to 


h{ 0 ,i,j,l) = 


rfj. 


-2e /-I 


2e 


dx[r {i X + (1 — x)j ) -I- 1] 2. 


(A33) 


Now, all we have to do is to take into account the combinations of the integrals just computed. 
The nontrivial two-loop diagram of the four-point function has the following result: 






27 




I2e2 


+ 2eC 


(^)i 


i,r 


(A34) 


The two-point vertex part two-loop diagram which interests us is 
'N + 2' 




[h{k, i, (T, 12 ) + 3h{k, i, 0, a, 12 ) + ?>Th{k, 0, i, a, p)]. (A35) 


We actually need the derivative of this diagram computed at zero external momenta. Let 
us hrst consider /3((A:, z, d,/i). Factorizing the mass, it can be written as follows 

d^^-^qihiqi k,ji +r,r) 




{ql + + 1) 


(A36) 


3 l=-oo 

The solution of the subdiagram is then substituted into this expression and yields 

d‘^~^qi 


h{k,i,a,i2) = /z 


1 


(1 


n CXD n 


Ji=-oo 




,2 ^2j2 ^ +r 2 (j^ + i)2 + 1] 

(A37) 
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where F^](qi + k,ji + i; r) was defined in Eq. fl90al) . The remaining of the compntation is 
very similar to what was worked ont for PBC and ABC] the reader is advised to consnlt 


Ref. ( 2^) for farther details. The derivative at zero external momentnm can be shown to 
be given by 


dh{kA]f) 


dk^ 




-2e , 


k^=0 86 


1 - - + eWo{r,r) 


where Wo{r,r) = Go{r,r) + Ho{r,r) — {r,r), 

“1 


Go{i;f) = “ 2 “^ 



dxdy{l — y)ln y{l — y)v + (1 — y)r ^ + 


yr 


-2 


0 JO 


a;(l — x) 


1 


Ho{r,r) = 2 



dxdy{l - y)fi {iy, \ y{l - y)F + {I - y)r 2 + 


yr 


-2 


(A38) 

,(A39a) 

(A39b) 


0 JO 


with FQ^\i]r) dehned in Eq. flhObll . We have to compnte k) 

scale ont the mass to write it as 

d^'-^qj^^q + k, i] f) 


a;(l — x)' 

. First, we 


fc2=0 


2-e 


h{k,t,r,d,y) = kfi J 
From the solntion of the snbdiagram we have the following intermediate result 


(A40) 


h{k,iJ;r) = 


k 


2-2e 


fll-- 


dx[x{l — x)] 2 


d^ ^qhiy + k, i] r) 


q2 + Cp + 1 


-lr~ 


_ 2 _ 

(g2 q- f'^p -)- 1) 

Let us define the integral 


*3 = 


d^ ^qhiq + k, i; f) 


qi _|_ y.2j2 _|_ Y 


[q + ky + CF + 1 


{q + ky + CF + 1 

(A41) 

(A42) 


Introduce a parameter y in order to write it in the form 


*3 = 



r(l + |) d^^ ^q dy y^ 


-1 


!+■ 


-. (A43) 


r(|) + 2kqy + k'^y + r2(j2(l - y) + Fy) + l-y + 

Next integrate over q and transform Sd-i into Sd using the prescription already stated. 
When we use the dehnition Eq. fl90cp together with our labor on F leads us to conclude 
that 


h{k,i,j-,r) = 


k 


2-2e 


rT(2 - f)r(-4 + e) 


1 - 


-1 


Py{l -y) + P{p{l -y) + Fy) + l-y + r) 


dx[x{l — a;)] 2 j dyy^ 


5-1 


. (A44) 


64 








































Taking the derivative at zero external momentum, performing the expansion in e, discarding 
0(e) contributions and using the dehnition Eq. fl90dp we get the desired expression, namely 


dl 3 {k,i,j;r) 


dk^ 


fc2=0 


= -h 


-2e 


r 

L4 


I dx J dy{l-y){r {j {l-y)+i y) + l-y + 




x(l — x) - 

(A45) 


Plugging the previous expression in the combinations appearing in the two-loop graph, we 
get to an expression which goes beyond the result from the simple PBC and ABC, or in 
other words 


dk-i (<o< 


fc2=0 


yi-‘^\N + 2) 
24e 

»1 /■! 


1 - - + eiyM(i,f) 
4 . 


W (z, r) = Wo(z; r) -|- 6r 
X (1 - J/) -y) + l-y + 


0 


dx I dy{l — y)(r‘^i‘^y + 1 — y + 

y 


(A46a) 


+ T dx dy 

Jo Jo 


x(l — x ). 


1 

\ — TT 



- 12 




x(l — x) 

(z, 0; r) + r (0, r)! • (A46b) 


The three-loop diagram for the two-point function has a similar systematics. We just need 

{k,i,0,d,ii) I 1 ^isik,0,^,a,^i) I 

\k^=0 &nQ \k^=0- 


A:2=0 


to compute a combination of the objects 
The hrst one can be checked to be given by the following expression: 

dh{k,i,d,fi) .1.. . qtWPz. 


dk^ 




e 

4 2 


(A47) 


where Jo ?/(l “+ (1 - 2/)^ ^ + ^ffz^ + 2 Jo dxdy(l - 

y)fi{iy, \Jy{^ — y)iJ + (1 ~ y)r~‘^ + — 4FQ(z;r). This is precisely equal to the con¬ 

tribution from PBC in Ref. 22|. We shall discuss the computation of the nondiagonal 
contribution. First, notice that I^ik, i,j, a, y) can be written in the form 

d‘^-^qll{q + k,j,f) 


h{k,i,j,a,y) = y V 


(A48) 


^2 I j^2.^2 I 

Replacing the value of the subdiagram already computed previously (with Sd already di¬ 
vided), we hnd 

d‘^~^q 


hik,i,j,a,y) = y 


2-3e 


r dx 


[g2 + r'^p + 1] (g + k)'^ + r'^f + 


-|-re 


d‘^ ^gF| (g + fc,j,r) 

g 2 _|_ ^ 2^2 _|_ I 


(A49) 
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Now utilize another Feynman parameter, integrate over q (transform the product 5'rf_ir(^^) 
into S'dr(|) as described in the text) and divide the result by Sd- Expand the remaining F 
functions in e and neglect 0{e^). Next, taking the derivative at zero external momenta we 
hnd 




dk^ 

y 


fc2=0 


= h 


-3e 


1 

Ye 




'F) 


0 


r / dy{l -y) dx r {i {l-y)+j y) + l-y 


(A50) 


x(l — x)- 

where j,r) was defined in Eq. (I90dp . The other terms involve the integration of 

products of the type I 2 I 2 and /2-f2- The former can be shown to be 0(e°) whereas the 
former are 0(e). Both terms which contribute a total of ten integrals in the computation of 
the three-loop graph of the two-point function are all regular and can be neglected in the 
evaluation of the singular part of the diagram. 

With the resources furnished in this Appendix, the reader is invited to check that the 
following expression holds 




fc2=0 


+ 2){N+ 8) 

162e2 


i-j + f f) 


(A51) 
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